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Abstract 

We demonstrate relationships between the classic Euclidean algorithm and many other fields 
of study, particularly in the context of music and distance geometry. Specifically, we show how 
the structure of the Euclidean algorithm defines a family of rhythms which encompass over forty 
timelines (ostinatos) from traditional world music. We prove that these Euclidean rhythms have 
the mathematical property that their onset patterns are distributed as evenly as possible: they 
maximize the sum of the Euclidean distances between all pairs of onsets, viewing onsets as 
points on a circle. Indeed, Euclidean rhythms are the unique rhythms that maximize this 
notion of evenness. We also show that essentially all Euclidean rhythms are deep: each distinct 
distance between onsets occurs with a unique multiplicity, and these multiplicies form an interval 
1, 2, . . . , A: — 1. Finally, we characterize all deep rhythms, showing that they form a subclass of 
generated rhythms, which in turn proves a useful property called shelling. All of our results for 
musical rhythms apply equally well to musical scales. In addition, many of the problems we 
explore are interesting in their own right as distance geometry problems on the circle; some of 
the same problems were explored by Erdos in the plane. 
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1 Introduction 



Polygons on a circular lattice, African bell rh ythms TouOSl ] . musical scales CEK99l |. sp allation 



neutron source accelerators in nuclear physics Bio03b| , linear sequences in mathematics LP92| 



mechanical words and str ingology in computer science dr awing digital straight lines in 

computer graphics calculating leap years in calendar design |HR04l . lAscOi |. and an ancient 

algorit hm for compu ting the greatest common divisor of two numbers, originally described by 
Euclid lEucSd iFraSGl ] — what do these disparate concepts all have in common? The short answer 
is, "patterns distributed as evenly as possible". For the long answer, please read on. 

Mathematics and music have been intimately intertwined since over 2,500 years ago when the 
famous Greek mathematician, Pythagoras of Samos (circa 500 B.C.), di scovered that the pleasing 



AshOal ]. Most of this in- 



experience of musical harmony is the result of ratios of small integers 
teraction between the two fields, however, has been in the domain of pitch and scales. For some 
histori cal sna pshots of this interaction, we refer the reader to H. S. M. Coxeter's delightful ac- 



Cox62l |. In mu sic theory, much attention has been devoted to the study of intervals used 



count ^ ^ 

in pitch scales For73l | . but relatively little work has been devoted to the analysis o f time d uration 
intervals of rhyt hm. Some notable recen t excep tions are the books by Simha Arom Aro9ll ]. Justin 
London jLon04| and Christopher Hasty Has97l | . 

In this paper, we study various mathematical properties of musical rhythms and scales that 
are all, at some level, connected to an algorithm of another famous ancient Greek mathematician, 
Euclid of Alexandria (circa 300 B.C.). We begin (in Section [2]) by showing several mathematical 
connections between musical rhythms and scales, the work of Euclid, and other areas of knowledge 
such as nuclear physics, calendar design, mathematical sequences, and computer science. In par- 
ticular, we define the notion of Euclidean rhythms, generated by an algorithm similar to Euclid's. 
Then, in the more technical part of the paper (Sections [3l l5.2p . we study two important properties 
of rhythms and scales, called evenness and deepness, and show how these properties relate to the 
work of Euclid. 

The Euclidean algorithm has been connected to music theory previously by Viggo Brun Bru64 1 . 
Brun used Euclidean algorithms to calculate the lengths of strings in musical instruments between 
two lengths I and 21, so that all pairs of adjacent strings have the same length ratios. In contrast, 
we relate the Euclidean algorithm to rhythms and scales in world music. 

Musical rhythms and scales can both be seen as two-way infinite binary sequences Tou02f | . In 
a rhythm, each bit represents one unit of time called a pulse (for example, the length of a sixteenth 
note), a one bit represents a played note or onset (for example, a sixteenth note), and a zero bit 
represents a silence (for example, a sixteenth rest). In a scale, each bit represents a pitch (spaced 
uniformly in log- frequency space) , and zero or one represents whether the pitch is absent or present 
in the scale. Here we suppose that all time intervals between onsets in a rhythm are multiples of a 
fixed time unit, and that all tone intervals between pitches in a scale are multiples of a fixed tonal 
unit (in logarithm of frequency). 

The time dimension of rhythms and the pitch dimension of scales have an intrinsically cyclic 
nature, cycling every measure and every octave, respectively. In this paper, we consider rhythms 
and scales that match this cyclic nature of the underlying space. In the case of rhythms, such 
cyclic rhythms are also called timelines, rhythmic phrases or patterns that are repeated throughout 
a piece; in the remainder of the paper, we use the term "rhythm" to mean "timeline" . The infinite 
bit sequence representation of a cyclic rhythm or scale is just a cyclic repetition of some n-bit string, 
corresponding to the timespan of a single measure or the log-frequency span of a single octave. To 
properly represent the cyclic nature of this string, we imagine assigning the bits to n points equally 
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spaced around a circle of circumference n McC9^. A rhythm or scale can therefore be represented 



as a subset of these n points. We use k to denote the size of this subset; that is, k is the number of 
onsets in a rhythm or pitches in a scale. For uniformity, the terminology in the remainder of this 
paper speaks primarily about rhythms, but the notions and results apply equally well to scales. 

In this paper, we use four representations of rhythms of timespan n. The first representation 
is the commonly used box-like representation, also known as the Time Unit Box System (TUBS), 
whic h is a se quence of n ' x 's and ' • 's where ' x ' represents an onset and ' • ' denotes a silence (a zero 
bit) [Toun2l |. This notation was used and taught in the West by Philip Harland at the University 
of California, Lo s Angel es, in 1962, and it was made popular in the field of ethnomusicology by 



James Koetting KoeTOj]. However, such box notation has been used in Korea for hundreds of 
years HK8lf |. The second representation of rhythms and scales we use is the clockwise distance 
sequence, which is a sequence of integers that sum up to n and represent the lengths of the intervals 
between consecutive pairs of onsets, measuring clockwise arc-lengths or distances around the circle 
of circumference n. The third representation of rhythms and scales writes the onsets as a subset of 
the set of all pulses, numbered 0, 1, . . . , n — 1, with a subscr ipt of n on the right-hand side of the 



subset to denote the timespan. Clough and Douthett jCD9ll | use this notation to represent scales. 
For example, the Cuban clave Son rhythm can be represented as[x - - x- - x- • - x- x- • •] 
in box-like notation, (3, 3, 4, 2, 4) in clockwise distance sequence notation, and |0, 3, 6, 10 , 12}i6 in 
subset notation. Finally, the fourth representation is a graphical clock diagram |Tou02[ | . such as 
Figure [H in which the zero label denotes the start of the rhythm and time flows in a clockwise 
direction. In such clock diagrams we usually connect adjacent onsets by line segments, forming a 
polygon. 

Even Rhythms. Consider the following three 12/8-time rhythms expressed in box-like notation: 
[x-x-x-x-x-x-], [x-x-xx-x-x-x], and [x • • • x x • • x x x • ]. It is intuitively clear 
that the first rhythm is more even (well spaced) than the second rhythm, and that the second 
rhythm is more even than the third rhythm. In fact, the second rhythm is the internationally most 
well known of all African timelines. It is traditi onally p layed on an iron bell, and is known on the 



world scene mainly by its Cuban name Bembe Tou03l | . Traditional rhythms tend to exhibit such 
properties of evenness to some degree. 

Why do many traditional rhythms display such evenness? Many are timelines (also sometimes 
called claves ) , that is, rhythms repeated throughout a piece that serve as a rhythmic reference 



point Uriod . Ort95 |. Often these claves have a call-and-response structure, meaning that the 



pattern is divided into two parts: the first poses a rhythmic question, usually by creating rhythmic 
tension, and the second part answers this question by releasing that tension. A good example of this 
structure is the popular clave Son [x - - x- - x- • - x- x- • •]. This clave creates such tension 
through syncopation, which can be found between the second and third onsets as well as between 
the third and fourth onsets. The latter is weak syncopation because the strong beat at position 8 lies 
half-way between the third and fourth onsets. (The strong beats of the underlying 4/4 meter (beat) 
occur at positions 0, 4, 8, and 12.) On the other hand, the former syncopation is strong bec ause the 
strong beat at position 4 is closer to the second onset than to the third onset Claves 
played with instruments that produce unsustained notes often use syncopation and accentuation 
to bring about rhythmic tension. Many clave rhythms create syncopation by evenly distributing 
onsets in contradiction with the pulses of the underlying meter. For example, in the clave Son, 
the first three onsets are equally spaced at the distance of three sixteenth pulses, which forms a 
contradiction because 3 does not divide 16. Then, the response of the clave answers with an offbeat 
onset, followed by an onset on the fourth strong beat of a 4/4 meter, releasing that rhythmic 
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(a) Shiko 



(b) Son 



(c) Soukous 




(d) Rumba (e) Bossa-Nova (f) Gahu 



Figure 1: The six fundamental African and Latin American rhythms which all have equal sum of pairwise 
geodesic distances; yet intuitively, the Bossa-Nova rhythm is more "even" than the rest. 



tension. 

On the other hand, a rhythm should not be too even, such as the example [x • x • x • x • x • 
X • ]. Indeed, in the most interesting rhythms with k onsets and timespan n, k and n are relatively 
prime (have no common divisor larger than 1). This property is natural because the rhythmic 
contradiction is easier to obtain if the onsets do not coincide with the strong beats of the meter. 
Also, we find that many claves have an onset on the last strong beat of the meter, as does the clave 
Son. This is a natural way to respond in the call-and-response structure. A different case is that 
of the Bossa-Nova clave [x - - x- - x- • - x- - x- •]. This clave tries to break the feeling of the 
pulse and, although it is very even, it produces a cycle that perceptually does not coincide with 
the beginning of the meter. 

This prevalence of evenness in world rhythms has led t o the study of rnathema tical measures of 
evenness in the new field of mathematical ethnomusicology Che02 . Tou04b . TouOSl ] , where they may 
help to identify, if not explain, cultural preferences of rhythms in traditional music. Furthermore, 
evenness in musical ch ords pla ys a significant role in the efficacy of voice leading as discussed in 
the work of Tymoczko TvmOd ] . 

The notion of maxim ally ev en sets with respect to scales repres ented o n a circle was introduced 
by Clough and Douthett CD9l[ | . According to Block and Douthett [BD94l | , Douthett and Entringer 
went furthe r by con structing several mathematical measures of the amount of evenness contained in 
a scale; see [BD94l . page 40]. One of their evenness measures simply sums the interval arc-lengths 
(geodesies along the circle) between all pairs of onsets (or more precisely, onset points). This 
measure differentiates between rhythms that differ widely from each other. For example, the two 
four-onset rhythms [x • • • x • • • x • • • x • • • ] and [x • x • x • • x ] yield evenness 



values of 32 and 23, respectively, reflecting clearly that the first rhythm is more evenly spaced than 
the second. Howeve r, the raeasure is too coarse to be useful for comparing rhythm timelines such 
as those stud ied in Tou02l . lTou03l | . For example, all six fundamental 4/4-time clave/bell patterns 
discussed in [Tou02l | and shown in Figure [T] have an equal pairwise sum of geodesic distances, 
namely 48, yet the Bossa-Nova clave is intuitively more even than, say, the Soukous and Rumba 
claves. 

Another distance measure that has been considered is the sum of pairwise chordal distances 
between adjacent onsets, measured by Euclidean distance between points on the circle. It can 
be shown that the rhythms maximizi ng this measure of evenness are precisely the rhythms with 
maximum possible area. Rappaport RapOSl ] shows that many of the most common chords and 
scales in Western harmony correspond to these maximum-area sets. This evenness measure is finer 
than the sum of pairwise arc- lengths, but it still does not distinguish half the rhythms in Figured) 
Specifically, the Son, Rumba, and Gahu claves have the same occurrences of arc-lengths between 
consecutive onsets, so they also have the same occurrences (and hence total) of distances between 
consecutive onsets. 

The measure of evenness we consider here is the sum of all pairwise Euclidean distances be- 

This measure is more 



BD9|]. 



tween points on the circle, as described by Block and Douthett 
discriminating than the others, and is therefore the preferred measure of evenness. For example, 
this measure distinguishes all of the six rhythms in Figure [H ranking the Bossa-Nova rhythm as 
the most even, followed by the Son, Rumba, Shiko, Gahu, and Soukous. Intuitively, the rhythms 
with a larger sum of pairwise chordal distances have more "well spaced" onsets. 

In Section IH we study the mathematical and computational aspects of rhythms that maximize 
evenness. We describe three algorithms that generate such rhythms, show that these algorithms 
are equivalent, and show that in fact the rhythm of maximum evenness is essentially unique. These 
results characterize rhythms with maximum evenness. One of the algorithms is the Euclidean- 
like algorithm from Section O proving that the rhythms of maximum evenness are precisely the 
Euclidean rhythms from that section. 



Deep Rhythms. Another important property of rhythms and scales that we study in this paper 
is deepness. Consider a rhythm with k onsets and timespan n, represented as a set of k points on a 
circle of circumference n. Now measure the arc- length/geodesic distances along the circle between 
all pairs of onsets. A musical scale or rhythm is Winograd-deep if every distance 1,2,..., [n/2j has 
a unique number of occurrences (called the multiplicity of the distance). For example, the rhythm 
[xxx • X •] is Winograd-deep because distance 1 appears twice, distance 2 appears thrice, and 
distance 3 appears once. 

The notion of deepness in scales w as introduced by Winograd in an oft-cited but unpublished 
class project rep ort from 1966 Wiri66l| . disseminated and further developed by the class instructor 
Gamer in 1967 [Gam67al . lGam67bl ]. and considered further in numerous papers and books, e.g., 
CEK99I . I Joh03l | . Equivalently, a scale is Winograd-deep if the number of onsets it has in common 



with ea ch of its cyclic shifts (rotations) is unique. This equ ivalence is the Common Tone Theorem 
Joh03, page 42], and it is originally described by Winograd Win66l | (who in fact uses this definition 
as his primary definition of "deep"). Deepn ess is o ne property of the ubiquitous Western diatonic 

and it captures some of the rich structure that 



,Tohn3l | 



12-tone major scale [x-x-xx-x-x- 
perhaps makes this scale so attractive. 

Winograd-deepness translates directly from scales to r hythrn s . For e xample, the diatonic major 
scale is equivalent to the famous Cuban rhythm Bembe PreSSl . Tou03 |. Figure [2] shows a graph- 
ical example of a Winograd-deep rhythm. However, the notion of Winograd-deepness is rather 
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Figure 2: A rhythm with k ~ 7 onsets and timespan n — 16 that is Winograd-deep and tlius Erdos-deep. 
Distances ordered by multiplicity from 1 to 6 are 2, 7, 4, 1, 6, and 5. The dotted line shows how the rhythm 
is generated by multiples of to = 5. 



restrictive for rhythms, because it requires half of the pulses in a timespan (rounded to a nearest 
integer) to be onsets. In contrast, for example, the popular Bossa-Nova rhythm [x • • x • • x • • • 
X • • X • • ] = {0, 3, 6, 10, 13}i6 pictures in Figure [1] has only five onsets in a timespan of sixteen. 
Nonetheless, if we focus on just the distances that appear at least once between two onsets, then 
the multiplicities of occurrence are all unique and form an interval starting at 1: distance 4 occurs 
once, distance 7 occurs twice, distance 6 occurs thrice, and distance 3 occurs four times. 

We therefore define a rhythm (or scale) to be Erdos-deep if it has k onsets and, for every 
multiplicity 1, 2, . . . , fc — 1, there is a nonzero arc- length/geodesic distance determined by t he point s 
on the circle with exactly that multiplicity. The same definition is made by Toussaint [Tou04a| | . 
Every Winograd-deep rhythm is also Erdos-deep, so this definition is strictly more general. 

To further clarify the difference between Winograd-deep and Erdos-deep rhythms, it is useful 
to consider which distances can appear. For a rhythm to be Winograd-deep, all the distances 
between 1 and k — 1 must appear a unique number of times. In contrast, to be an Erdos-deep 
rhythm, it is only required that any distance that appears must have a unique multiplicity. Thus, 
the Bossa-Nova rhythm is not Winograd-deep because distances 1,2 and 5 do not appear. 

The property of Erdos deepness involves only the distances between points in a set, and is 
thus a feature of distance geometry — in this_case, in the discrete space of n points equally spaced 
around a circle. In 1989, Paul Erdos jErdSgl ] considered the analogous question in the plane, 
asking whether there exist n points in the plane (no three on a line and no four on a circle) such 
that, for every i = 1, 2, . . . , n — 1, there is a distance determined by these points that occurs exactly 
i times. Solut ions have been found for n between 2 and 8, but in general the problem remains open. 
Palasti PalBDi] considered a variant of this problem with further restrictions — no three points form 
a regular triangle, and no one is equidistant from three others — and solved it for n = 6. 

In Section m we characterize all rhythms that are Erdos-deep. In particular, we prove that all 
deep rhythms, besides one exception, are generated, meaning that the rhythm can be represented 
as {0, m, 2m, . . . ,(k — l)m}n for some integer m, where all arithm etic is m odulo n. In the context 
of s cales, th e concept of "generated" was defined by Wooldridge Woo93( | and used by Clough et 
al. [CEK99 |. For example, the rhythm in Figure [2] is generated with m = 5. Our ch aracteriz ation 
generalizes a similar characteriz ation fo r Winograd-deep scales proved by Winograd Win6d | . and 
independently by Clough et al. [CEK99 |. 

In the pitch domain, generated scales are very common. The Pythagorean tuning is a good 
example: all its pitches are generated from the fifth of ratio 3 : 2 modulo the octave. Anothe r 
example is the equal-tempered scale, which is generated with a half-tone of ratio \^ }Bar04a |. 
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Generated scales are also of interest in the theory of the well- formed scales |Car98l ] . 

Generated rhythms have an interesting property called shellability. If we remove the "last" 
generated onset 14 from the rhythm in Figure [21 the resulting rhythm is still generated, and this 
process can be repeated until we run out of onsets. In general, every generated rhythm has a 
shelling in the sense that it is always possible to remove a particular onset and obtain another 
generated rhythm. 

Shellings of rhythms play an important role in musical improvisation. For example, most African 
drumming music consists of rhythms operating on three different strata: the unvarying timeline 
usually provided by one or more bells, one or more rhythmic motifs played on drums, and an 
improvised solo (played by the lead drummer) riding on the other rhythmic structures. Shellings 
of rhythms are relevant to the improvisation of solo drumming in the context of such a rhythmic 
background. The solo improvisation must respect the style and feeling of the piece which is usually 
determined by the timeline. A common technique to achieve this effect is to "borrow" notes from 
the timeline, and to alternate between pl aying subsets of notes from the timeline and fro m othe r 
rhythms that interlock with the timeline AnkQTI . In the words of Kofi Agawu [AgaSfil ]. 

"It takes a fair amount of expertise to create an effective improvisation that is at the same time 
stylistically coherent". The borrowing of notes from the timeline may be regarded as a fulfillment 
of the requirements of style coherence. 

Of course, some subsets of notes of a rhythm may be better choices than others. For example, 
it seems reasonable that, if a rhythm is deep, one should select subsets of the rhythm that are also 
deep. Furthermore, a shelling seems a natural way to decrease or increase the density of the notes 
in an improvisation that respects these constraints. For example, in the Bembe bell timeline [x • 
X • X X • X • X • x], which is deep, one possible shelhng is [x • x • x x • x • x • • ], [x • x • x • • 
x-x - •], [x-x - • • - x- x- •], [x-x - • • - x- • • •]. All five rhythms sound good and are 
stylistically coherent. To our knowledge, shellings have not been studied from the musicological 
point of view. However, they may be useful both for theoretical analysis as well as providing formal 
rules for "improvisation" techniques. 

One of the consequences of our characterization that we obtain in Section [5] is that every Erdos- 
deep rhythm has a shelling. More precisely, it is always possible to remove a particular onset that 
preserves the Erdos-deepness property. Finally, to tie everything together, we show that essentially 
all Euclidean rhythms (or equivalently, rhythms that maximize evenness) are Erdos-deep. 



2 Euclid and Evenness in Various Disciplines 

In this section, we first describe Euclid's classic algorithm for computing the greatest common 
divisor of two integers. Then, through an unexpected connection to timing systems in neutron 
accelerators, we see how the same type of algorithm can be used as an approach to maximizing 
"evenness" in a binary string with a specified number of zeroes and ones. This algorithm defines an 
important family of rhythms, called Euclidean rhythms, which we show appear throughout world 
music. Finally, we see how similar ideas have been used in algorithms for drawing digital straight 
lines and in combinatorial strings called Euclidean strings. 



2.1 The Euclidean Algorithm for Greatest Common Divisors 

The Euclidean algorithm for computing the greatest common divisor of two integers is one of the 
oldest known algorithm s (circa 300 B .C.). It was first describe d by E uclid in Proposition 2 of 
Book VII of Elements [ EucSS . Fra56 |. Indeed, Donald Knuth Knu98l | calls this algorithm the 
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"granddaddy of all algorithms, because it is the oldest nontrivial algorithm that has survived to 
the present day" . 

The idea of the algorithm is simple: repeatedly replace the larger of the two numbers by their 
difference until both are equal. This final number is then the greatest common divisor. For example, 
consider the numbers 5 and 13. First, 13 — 5 = 8; then 8 — 5 = 3; next 5 — 3 = 2; then 3 — 2 = 1; 
and finally 2 — 1 = 1. Therefore, the greatest common divisor of 5 and 13 is 1; in other words, 5 
and 13 are relatively prime. 



The algorithm can also be described succinctly in a recursive manner as follows CLRSOll ] . Let 
k and n be the input integers with k < n. 



Algorithm EuCLlD(fc,n) 

1. if = then return n 

2. else return EuCLiD(n mod k, k) 



Running this algorithm with k = 5 and n = 13, we obtain Euclid(5, 13) = Euclid(3, 5) = 
Euclid(2,3) = Euclid(1,2) = Euclid(0, 1) = 1. Note that this division version of Euclid's 
algorithm skips one of the steps (5, 8) made by the original subtraction version. 

2.2 Evenness and Timing Systems in Neutron Accelerators 

One of our main musical motivatations is to find rhythms wi th a specified timespan and number 



of onsets that maximize evenness. Bjorklund [Bio03bl . lBio03a[ | was faced with a similar problem of 



maximizing evenness, but in a different context: the operation of components such as high-voltage 
power supplies of spallation neutron source (SNS) accelerators used in nuclear physics. In this 
setting, a timing system controls a collection of gates over a time window divided into n equal- 
length intervals. (In the case of SNS, each interval is 10 seconds.) The timing system can send 
signals to enable a gate during any desired subset of the n intervals. For a given number n of time 
intervals, and another given number A; < n of signals , the pr oblem is to distribute the pulses as 



evenly as possible among these n intervals. Bjorklund [Bio03bl | represents this problem as a binary 
sequence of k ones and n — k zeroes, where each bit represents a time interval and the ones represent 
the times at which the timing system sends a signal. The problem then reduces to the following: 
construct a binary sequence of n bits with k ones such that the k ones are distributed as evenly as 
possible among the (n — k) zeroes. 

One simple case is when k evenly divides n (without remainder), in which case we should place 
ones every n/k bits. For example, if n = 16 and A; = 4, then the solution is [1000100010001000]. 
This case corresponds to n and k having a common divisor of k. More generally, if the greatest 
common divisor between n and k is g, then we would expect the solution to decompose into g 
repetitions of a sequence of n/g bits. Intuitively, a string of maximum evenness should have this 
kind of symmetry, in which it decomposes into more than one repetition, whenever such symmetry is 
possible. This connection to greatest common divisors suggests that a rhythm of maximum evenness 
might be computed using an algorithm like Euclid's. Indeed, Bjorklund's algorithm closely mimics 
the structure Euclid's algorithm, although this connection has never been mentioned before. 

We describe Bjorklund's algorithm by using one of his examples. Consider a sequence with 
n = 13 and k = 5. Because 13 — 5 = 8, we start by considering a sequence consisting of 5 ones 
followed by 8 zeroes which should be thought of as 13 sequences of one bit each: 

[1][1][1][1][1][0][0][0][0][0][0][0][0] 
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If there is more than one zero the algorithm moves zeroes in stages. We begin by taking zeroes 
one at a time (from right to left), placing a zero after each one (from left to right), to produce five 
sequences of two bits each, with three zeroes remaining: 

[10] [10] [10] [10] [10] [0] [0] [0] 

Next we distribute the three remaining zeros in a similar manner, by placing a [0] sequence after 
each [10] sequence: 

[100] [100] [100] [10] [10] 

Now we have three sequences of three bits each, and a remainder of two sequences of two bits each. 
Therefore we continue in the same manner, by placing a [10] sequence after each [100] sequence: 

[10010] [10010] [100] 

The process stops when the remainder consists of only one sequence (in this case the sequence 
[100]), or we run out of zeroes (there is no remainder). The final sequence is thus the concatenation 
of [10010], [10010], and [100]: 

[1001010010100] 

We could proceed further in this process by inserting [100] into [10010] [10010]. However, 
Bjorklund argues that, because the sequence is cyclic, it does not matter (hence his stopping rule). 
For the same reason, if the initial sequence has a group of ones followed by only one zero, the 
zero is consi dered a s a remainder consisting of one sequence of one bit, and hence nothing is done. 



Bjorklund [Bjo03bl ] shows that the final sequence may be computed from the initial sequence using 
0{n) arithmetic operations in the worst case. 

A more convenient and visually appealing way to implement this algorithm by hand is to perform 
the sequence of insertions in a vertical manner as follows. First take five zeroes from the right and 
place them under the five ones on the left: 

111110 


Then move the three remaining zeroes in a similar manner: 

11111 



Next place the two remainder columns on the right under the two leftmost columns: 

1 1 1 




1 1 


Here the process stops because the remainder consists of only one column. The final sequence is 
obtained by concatenating the three columns from left to right: 

1001010010100 

Bjorklund's algorithm applied to a string of n bits consisting of k ones and n — k zeros has 
the same structure as running Euclid(A;, n). Indeed, Bjorklund 's algor ithm uses the repeated 
subtraction form of division, just as Euclid did in his Elements EucSd ]. It is also well known 



that applying the algorithm Euclid(A;, n) to two 0{n) bit number s (binary sequences of length n 



causes it to perform 0{n) arithmetic operations in the worst case CLRSOl ]. 
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Figure 3: (a) The Euclidean rhythm E{3, 8) is the Cuban tresillo. (b) The Euchdean rhythm £'(5, 8) is the 
Cuban cinquillo. 



2.3 Euclidean Rhythms 



The binary sequences generated by Bjorklund's algorithm, as described in the preceding, may 
be considered as one family of rhythms. Furthermore, because Bjorklund's algorithm is a way 
of visualizing the repeated-subtraction version of the Euclidean algorithm, we call these rhythms 
Euclidean rhythms. We denote the Euclidean rhythm by E{k,n), where k is the number of ones 
(onsets) and n (the number of pulses) is the length of the sequence (zeroes plus ones). For ex- 
ample, £'(5,13) = [1001010010100]. The zero-one notation is not ideal for representing binary 
rhythms because it is difficult to visualize the locations of the onsets as well as the duration of 
the inter-onset intervals. In the more iconic box notation, the preceding rhythm is written as 
£^(5,13) = [x • • X • X • • X • X • •]. 



The rhythm E{5, 13) is in fact used in Macedonian music [AroOj], but having a timespan of 13 
(and defining a measure of length 13), it is rarely found in world music. For contrast, let us consider 
two widely used values of k and n; in particular, what is £(3,8)? Applying Bjorklund's algorithm 
to the corresponding sequence [11100000], the reader may easily verify that the resulting Euclidean 
rhythm is £(3,8) = [x • • x • • x • ]. Figure [3ja) shows a clock diagram of this rhythm, where the 
numbers by the sides of the triangle indicate the arc-lengths between those onsets. 

The Euclidean rhythm £(3, 8) is one of the most famous on the planet. In Cuba, it goes by 
the name of the tresillo, and in the USA, it is often called the Habanera rhythm. It was used in 
hundreds of rockabilly songs during the 1950's. It can often be heard in earl y rock - and-roll hits i n 
the left-hand patterns of the piano, or played on the string bass or sa xophon e Bre99l . Flo99l . Mor96 ]. 
A good example is the bass rhythm in Elvis Presley's Hound Dog Bre99|. The tresillo pattern is 
also found widely in West African tra ditional music. For example, it is played on the atoke bell in 
the Sohu, an Ewe dance from Ghana KauBOl ] . The tresillo can also be recognized as the first bar 
(first eight pulses) of the ubiquitous two-bar clave Son shown in Figure [U^b) . 

In the two examples £(5, 13) and £(3,8), there are fewer ones than zeros. If instead there are 
more ones than zeros, Bjorklund's algorithm yields the following steps with, for example, k = 5 
and n = 8: 



[1 1 1 1 1 0] 
[10] [10] [10] [1] [1] 
[101] [101] [10] 
[10 110 110] 
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Figure 4: These two rhythms are instances of the same rhythm necklace. 



The resulting Euclidean rhythm is E{5, 8) = [x • xx • xx •]. Figure[3)[|b) shows a clock diagram 
for this rhythm. It is another famous rhythm on the w orld sc ene. In Cuba, it goes by the name of 
the cinquillo and i t is inti mately related to the tresillo Flo99f | . It has been used in jazz throughout 



the 20th century |R,ah9fil ]. and in roc kabilly music. For example, it is the hand-clapping pattern 



in Elvis Presley's Hound Dog BreOOl ] . The cinquillo pattern is also widely used in West African 



traditional music jRahSTI . lTou02l |. as well as Egyptian [Hag03] and Korean [HK8l[ | music. 

We show in this paper that Euclidean rhythms have two important properties: they maximize 
evenness and they are deep. The evenness property should come as no surprise, given how we 
designed the family of rhythms. To give some feeling for the deepness property, we consider the 
two examples in Figure [3l which have been labeled with the distances between all pairs of onsets, 
measured as arc-lengths. The tresillo in Figure [3]^a) has one occurrence of distance 2 and two 
occurrences of distance 3. The cinquillo in Figure [3l^b) contains one occurrence of distance 4, two 
occurrences of distance 1, three occurrences of distance 2, and four occurrences of distance 3. Thus, 
every distance has a unique multiplicity, making these rhythms Erdos-deep. 

2.4 Euclidean Rhythms in Traditional World Music 

In this section, we list all the Euclidean rhythms found in world music that we have collected 
so far, restricting attention to those in which k and n are relatively prime. In some cases, the 
Euclidean rhythm is a rotated version of a commonly used rhythm. If a rhythm is a rotated version 
of another, we say that they are instances of the same necklace. Thus a rhythm necklace is a 
clockwise distance sequence that disregards the starting point in the cycle. Figure [4] illustrates an 
example of two rhythms that are instances of the same necklace. 

Rhythms in which k and n have a common divisor larger than 1 are common all over the 
planet in traditional, classical, and popular genres of music. For example, i?(4, 12) = [x • • x • • 
X • • X • • ] is the 12/8-time Fandango clapping pattern in the Flameii co music of southern Spain, 



where 'x' denotes a loud clap and ' • ' denotes a soft clap DBFG"'"04| ]. However, the string itself 
is periodic: -E(4, 12) has period 3, even though it appears in a timespan of 12. For this reason, we 
restrict ourselves to the more interesting Euclidean rhythms that do not decompose into repetitions 
of shorter Euclidean rhythms. We are also not concerned with rhythms that have only one onset 
([x • ], [x • • ], etc.), and similarly with any repetitions of these rhythms (for example, [x • x • ]). 

There are surprisingly many Euclidean rhythms with k and n relatively prime that are found 
in world music. The following list includes more than 40 such rhythms uncovered so far. 

£'(2,3) = [x X • ] = (12) is a common Afro-Cuban drum pattern when started on the second 
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onset as in [x • x]. For example, it is the conga rhythm of the (6/8)-tii ne Swin g Tunibao K1697l | . 
It is c ommon in Latin American music, as for example in the Cueca vdL95l |. and the cows de 
clave RodOTIl . It is common in Arabic music, as for example in the Al Tder rhythm of Nu- 



bia iHagn.'^. It is also a rhythmic pattern of the Drum Dance of the Slavey Indians of Northern 



Canada 



Asc75l |. 



Africa AroO 



X • X 



] = (23) is a rhythm found in Greece, Namibia, Rwan da and Central 
It is also a 13th century Persian rh ythm c alled Khafif-e-ramal Wri78l | . as well 



as the rhythm of the Macedonian dance Makedonka ISin74| . Tchaikovsky used it as the metric 
pattern in the second movement of his Symphony No. 6 [ Kei91 |. Started on the second on set as i n 
[x • • X • ] it is a rhythm found in Central Africa, Bulgaria, Turkey, Turkestan and Norway [AroOj]. 
It is also the metric pattern of Dave Brubeck's Take Five, as well as Mars from The Pla nets b y 
Gustav Hoist [Kei9l|. Both starting points determine metric patterns used in Kor ean mu sic [HK8l[ |. 

£'(3,4) = [x X X • ] = (112) is a pattern used in the Baiao rhythm of Brazil Uri93l | . as well as 
the polos of Bali MonSSi ] . Started on the second onset, as in [x x • x], it is the Catarete rhythm 
of the indigenous people of Brazil |Uri93l | . Started on the third onset, as in [x • xx], it is the 
archetypa l patte rn of the Cumbia from Colombia [Man85| | . as well as a Calyps o rhyt hm from 
Trinidad Eva66 |. It is also a 13th century Persian rhythm c alled K halif-e-saqhil WriTSi ] . as well 
as the trochoid choreic rhythmic pattern of ancient Greece jMat85l | . Started on the silent note 
(anacrusis), as in [• x x x ], it is a popular flamenco rhythm used in the Taranto, the Tiento, the 
Tango, and the Tanguillo Gam02l | . It is also the R umba c lapping pattern in flamenco, as well as a 
second pattern used in the Baiao rhythm of Brazil Uri93l | . 

£'(3, 5) = [x • X • x] = (221), when s tarted on the second onset, is another 13th century Persian 



rhythm by the name of Khafif-e-ramal Wri78 |. as well as a Ro manian folk-dance rhythm PC69| 



and the Sangsa Pydlgok drum pattern in Korean music [HK8ll |. 

E{2>, 7) = [x • X • X • • ] = (223) is a rhythm found in Greece, Tur kestan, Bulgaria, and North- 
ern Sudan [Aro04 ] . It is the Ddwer turan rhyth mic pattern of Turkey Hag03| | . It is the Ruchenitza 



rhythm used i n a B ulgarian folk-dance 
Eleno Mome 



Pre83l |. 



Sin7|. 



as well as the rhythm of the Macedonian dance 
It is also the rhythmic pattern of Dave Brubeck's Unsquare Dance, and 
Pink Fl oyd's M oney Kei91[. Started on the second onset, as in [x • x • • x • ], it is a Serbian 
rhythm AroOj]. St arted o n the third onset, as in [x • • x • x • ], it is a rhyt hmic pa ttern found in 
Greece and Turkey Aro04l | . In Yemen it goes by the n ame of Daasa al zreir [Hag03'l. It is also the 
rhythm of th e Maced onian dance Tropnalo Oro [SinTj], the rhythm for the Bulgarian Makedonsko 
Horo dance Wad04f |. as well as the meter and clapping pattern of the ttvra tdl of North Indian 
music ClaOol ^ 

E{3, 8) = [x • • X • • X • ] = (332) is the Cuban tresi llo patt ern discussed in the preceding j Flo99l | . 
the most important tradi tional b luegrass banjo rhythm Kei9l[ |, as wel l as the Mai metal-blade pat- 
tern of the Aka Pygmies Aro9ll |. It is com mon in . West Africa Bra59[l and many other parts o f the 
world such as Greece and Northern Sudan Aro04 1 . Curt Sachs Sac53l | and Willi Apel Ape6Cll | con- 
sider it to be one of the most important rhythms in Renaiss ance music. Indeed, it dates back to the 
Ancient Greeks wh o called it the dochmiac pattern Bra59l | . In India it is one of the 35 siiladi talas 
of Karnatak music LonOj]. Started on the second onset, it is a drum pattern of the Samhy on TodUri 
Korean instrumental music HK8l[ | . It is also found in Bulgaria and Turkey Aro9l[ | . Start ed on th e 
third onset, it is the Nandon Bawaa bell pat tern of t he Dagarti people of northwest Ghana HAD95l |. 
and is also found in Namibia and Bulgaria [Aro9ll ]. 

£'( 3, 11) = [x • • • X • • • X • • ] = (443) is the metric pattern of the savdri tdl of North Indian 
music |ClaOOl |. 

£(3,14) = [x • • • • X • • • • X • • •] = (554) is the clapping pattern of the rf/iamar iaZ of North 
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Indian music 
i?(4,5) = 



ClaOd |. 



[X X X X - J = (1112) is the rhythmic pattern of the Mirena rhyt hm of G reece Hag03| | . 
Started on the fourth onset, as in [x • x x x], it is the Tik rhythm of Greece Hag03l | . 

£'(4,7) = [x • X • X • x] = (2221) is another Ruchenitza Bulgarian folk -dance rhythm |Pre83l | . 
Started on the third onset, it is th e Kala mdtianos Greek dance rhythm jHag03l ] . as well as the 
Shaigie rhythmic pattern of Nubia jHag03 |. Star ted on the fourth (last) onset, it is the rhythmic 



Hag03l |. 



pattern of the Dar daasa al mutawasit of Yemen 

£(4,9) = [x-x-x-x - •] = (2223) is the Aksak rhythm of Turkey [Bra5l|l (also found in 
Greece) as well as the rhythi n of th e Macedonian dance Kambani Bijat Oro Sin74| and the Bul- 



garian dance Daichovo Horo Ric04l | . In Bulgarian music fast tunes with this metric patte rn are 
called Dajchovata whereas slow tunes with this same pattern are called Samokovskata fRicSOl] . It is 
the rhythmic ostinato of a lullaby discovered by Simha Arom in south-western ZaireJAroM. It is 
the metric pattern used by Dave Brubeck in his well known piece Rondo a la Turk [Kei9l|. Whe n 
it is started on the second onset, as in [x • x • x • • x • ], it is found in Bulgaria and Serbia fArpOi]. 
Started on the third onset, as in [x • x • • x • x • ], it is f ound in Bulgaria and Greece [Aro04l ] . 



It is the rhythm of the Macedonian dance Devojce SinTJ]. Finally , when started on the fourth 
onset, as in [x • • x • x • x • ], it is a rhythm found in Turkey AroOj], and is the metric pattern of 
Strawberry Soup by Don Ellis [Kei9l[. 

£( 4, 11) = [x - - x- - x- •x-] = (3332) is the metric patt ern of t he Dhruva tala of Southern 



Kei9l[ | . When it is started on the 



India LonOJ]. It is also used by Frank Zappa in Outside Now 
third onset, as in [x • • x • x • • x • • ], it is a Serbian rhythmic pattern Aro04l. Whe n it is started 
on the fourth (last) onset it is the Daasa al 
£(4,15) = [x • • • X • • • X • • • X • • 1 = 



biri rhythmic pattern of Yemen |llag03l ] . 
(4443) is the metric pattern of the pancam savdri tal 
of North Indian music 

£(5, 6) = [x X X X X • ] = (11112) yields the York-Samai pattern, a pop ular Ara bic rhythm Sta88l |. 
It is also a handclapping rhythm used in the Al Medemi songs of Oman [EMF9nl |. 

£( 5, 7) = [x • XX • xx] = (21211) is the Naw akhat p attern, another popular Arabic rhythm 
Stasi]. In Nubia it is called the Al Noht rhythm |llag03l |. 

£(5, 8) = [x • XX • XX • ] = (21212) is the C uban ci nquillo pattern discussed in the preced- 
ing [f1o99|, the Malfuf rhythmic pattern of Egypt Hag03l | . as well as the Korean Nong P'yon d[ rum 
pattern [HK8l[ |. Started on the second onset, it is a popular Middle Eastern rhythm WadOj], as 



Che79t . the Spanish 



well as the Ti mini rhythm of Senega l, the A dzoqbo dance rhythm of Benin ^ 

Tanao EvaGSl] . the Maksum of Egypt HagOS], and a 13th century Persian rhythm, the Al-s aqhil-a l- 
' Wri78| . When it is started on the third onset it is the Miisemmen rhythm of Turkey Bek05 |. 



sam 



When it is started on the fourth onset it is the Kromanti rhythm of Surinam. 

£(5,9) = [x • X • X • X • x] = (22221) is a popular Arabic rhythm called Agsag-Saniai 



Sta88|]. 



Started on the second onset, it is a drur n patt ern used by the Venda in South Africa Rah87l |. as 
well as a Rumani an folk- dance rhythm PC69l |. It is also the rhyth mic pat tern of the Sigaktistos 
rhythm of Greece Hag03], and the Samai aktsak rhythm of Turkey flIaeo3]. Started on the third 
onset, it is the rhythmic pattern of the Nawahiid rhythm of Turkey [Hag03]. 

£(5, 11) = [x • X • X • X • X • • ] = (22223) is the metric pattern of the SavarT tala use d in th e 
Hindustani music of India Lon04l | . It is also a rhythmic pattern used in Bulgaria and Serbia Aro04i | . 
In Bulgaria is is used in the Kopan itsa |Ric04 1 . This metric pattern has been used by Moussorgsky 
in Pictures at an Exhibit ion | Kei9ll| . Started on the third onset, it is the rh ythm of the Macedonian 
dance Kalajdzijsko Oro Sin74| - and it appears in Bulgarian music as well Aro04 |. 

£(5,12) = [x • • X • X • • X • X • ] = (3232 2 ) is a c ommon rhythm played in the Central African 



Republic by the Aka Pygmies Aro9ll . IChe02l . |CT03| . It is also the Venda clapping pattern of a 
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South African children's song Pre83l | . and a rhythm pattern used in Macedonia Aro04 ] . Starte d 
on the second onset, it is the Columbia bell pattern popular in Cuba a,nd Wes t Africa K1697l | . 
as well as a d rummi ng pattern used in the Chakacha dance of Kenya jBar04bl | . and also used 
in Macedor iia Aro04l] . Started on the third onset, it is the Bemba bell patterr i used in Northern 
Zimbabwe Pre83l |. and the rhythm of the Macedonian dance Ibraim Odza Pro Sin74 ] . Started on 
the fourth onset, it is the Fum e Fume bell pattern popular in West Africa [K1697f 7 and is a rhythm 
used in the former Yugoslavia AroOj]. Finally, when started on the fifth onset it is the Salve bell 
patter n used in the Dominican Republic in a rhythm called Canto de Vela in ho nor of t he Virgin 
Mary Far92l | . as well as the drum rhythmic pattern of the Moroccan Al Kuddm [HagOSl ]. 

E{5, 13) = [x - - x- x- - x- x- •] = (32323) is a Macedonian rhythin which is also played by 
starting it on the fourth onset as follows: 

£■(5, 16) = [x • • X • • X • • X • • X • • 



AroOj]. 



x-x- -X- -x-x- 

- ] = (33334) is the Bossa-Nova rhythm necklace of Brazil. 
The ac tual Bo ssa-Nova rhythm usually starts on the third onset as follows: [x - - x- - x- - - x- - 
X - - ] Tou02| ■ However, other st arting places are also documented in world music practices, such 
as[x - - X- - X- - X- - - X- - ] [Beh73^ . 

£'(6, 7) = [x X X X X X - ] = (111112) is the Pontakos rhythm of Greece when started on the sixth 
(last) onset |Hag03l | . 

E{6, 13 ) = [x - x- x- x- x- x- - ] = (222223) is the rhythm of the Macedonian dance Mama 
Con e pita Sin74l ]. Started on the third onset, it is the rhythm of the Macedonian dance Postupano 
Oro [Sin74| , as well as the Krivo Plovdivsko Horo of Bulgaria Ric04l | . 

£(7,8) = [xxxxxxx •] = (1111112), when started on the seventh (last) onset, is a typical 
rhythm played on the Bendir (frame drum) , and used in the accompaniment of songs of the Tuareg 
people of Libya [Sta88|]. 

E{7, 9) = [x - X X X - X X x] = (2112111) is the Bazaragana rhythmic pattern of Greece tHag03l]. 
£(7,10) = [x-xx-xx-xx] = (2121211) is the Len/;; /aMie rhythmic pattern of Turkey [Hag03(]. 
£(7,12) = [x-xx-x-xx-x-] = (2122122) is a common West Af rican b ell pattern. For 
example, it is used in the Mpre rhythm of the Ashanti people of Ghana Tou03|. Started on the 
seventh (last) onset, it is a Yoruba be ll patt ern of Nigeria, a Babenzele pattern of Central Africa, 
and a Mende pattern of Sierra Leone [Sto05l |. 

£(7, 15) = [x - x-x-x-x-x-x - - ] = (2222223) is a Bulgarian rhythm when started on the 



third onset Aro04 1 . 

^(7,16) = [x - - X - X - X 
The actual Samba rhythm is 



- X - X 
"x - X - 



■ X - ] = (3223222) is a Samba rhythm necklace from Brazil. 
X - X - X - - X - X - ] ob tained by starting £(7, 16) on the 



last onset, and it coincides with a Macedonian rhythm AroOJ]. When £(7, 16) is started on the 
fifth onset it is a clapping pattern from Ghana Pre83|. Whe n it is started on the second onset it 
is a rhythmic pattern found in the former Yugoslavia jAroo3 |. 

= (3232322) is a Macedonian rhythm when started 



X - X 



X - X 



X - X 



Sin7|. 



x-x - - X - X 
X - X - X - X - 



£(7,17) = [x 
on the second onset 

£(7,18) = [x - - X - X - 
£(8, 17) = [x - X - X - X 
is also started on the fifth onset Aro04i | . 

£(8, 19) = [x - - X - X - X - - X - X - X 
when started on the second onset 
£(9,14) = fx-xx-xx-xx-xx 



= (3232323) is a Bulgarian rhythmic pattern jAro04 1 . 
(22222223) is a Bulgarian rhythmic pattern which 

] = (32232232) is a Bulgarian rhythmic pattern 



AroOj]. 



^ = (212121 212), w hen started on the second onset, is the 
rhythmic pattern of the Tsofyan rhythm of Algeria Hag03l | . 

£(9, 16) = [x - X X - X - X - X X - X - X - ] = (212221222) is a rhythm necklace used in the Cen- 
tral African Republic Aro9l[. When it is started on the second onset it is a bell pattern of the 
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Luba people of Con go Mus02ll . When it is started on the fourth onset it is a rhyt hm pl ayed in West 
and Central Africa jFlo99l |. as well as a cow-bell pattern in the Brazilian samba [Sol96l | . When it is 
started on the pe nultima te onset it is the bell pattern of the Ngbaka-Maibo rhythms of the Central 
African Republic [Aro9ll |. 

E{9, 22) =[x • - x- x- - x- x- • X • X • • x • x • ] = (323232322) is a Bulgarian rhythmic pat- 



tern when started on the second onset AroOJ]. 

£■(9,23) =[x • - x- x- - x- x- - x- x- - x- x- •] = (323232323) is a Bulgarian rhythm [Aro04l | 
£'(11,12) =[x X X X X X X X X X X •] = (11111111112), when started on the second onset, is the 

drum pattern of the Rahmdni (a cylind rical double-headed drum) used in the Sot silam dance from 

Mirbdt in the South of Oman |EMF9d |. 

£(11, 24) =[x • - x- x- x- x- x- -x-x • x • x • x • ] = (32222322222) is a rhythm necklace 



of the Aka Pygmies of Central Africa Aro9l[ |. It is u sually started on the seventh onset. Started 
on the second onset, it is a Bulgarian rhythm Aro04 1 . 

£(13, 24) =[x - xx-x-x-x-x-xx-x-x • x • x • ] = (2122222122222) is another rhythm 



necklace of the Aka Pygmies of the upper Sangha Aro9ll ]. Started on the penultimate onset, it is 
the Bobangi metal-blade pattern used by the Aka Pygmies. 

£(15, 34) = [x- - x- x- x- x- - x- x- x- x- - x- x- x- x- - x • x • ] = (322232223222322) 



is a Bulgarian rhythmic pattern when started on the penultimate onset AroOJ]. 



2.5 Aksak Rhythms 

Euclidean rhythms are closely related to a family of rhythms knowi i as aksak rhythms, which 
have be en stud ied from the combinatori al poin t of view for some time Bra5ll . ICle94l . lAro04l | . Bela 
Bartok BarSlI ] and Constantin Brailoiu Bra5ll |. respectively, have used the terms Bulgarian rhythm 
and aksak to refer to those meters that use units of durations 2 and 3, and no other durations. 
Furthermore, the rhythm or me ter mu st contain at least one duration of length 2 and at least 
one duration of length 3. Arom [Aro04l | refers to these durations as binary cells and ternary cells, 
respectively; 



Arom Aro04i | generated an inventory of all the theoretically possible aksak rhythms for values 
of n ranging from 5 to 29, as well as a list of those that are actually used in traditional world 
music. He also proposed a classification of these rhythms into several classes, based on structural 
and numeric properties. Three of his classes are considered here: 

1. An aksak rhythm is authentic if n is a prime number. 

2. An aksak rhythm is quasi-aksak if n is an odd number that is not prime. 

3. An aksak rhythm is pseudo-aksak if n is an even number. 

A quick perusal of the Euclidean rhythms listed in the preceding reveals that aksak rhythms are 
well represented. Indeed, all three of Arom's classes (authentic, quasi-aksak, and pseudo-aksak) 
make their appearance. There is a simple characterization of those Euclidean rhythms that are 
aksak. ^From the iterative subtraction algorithm of Bjorklund it follows that \in = 2k all cells are 
binary (duration 2). Similarly, iin = 'ik all cells are ternary (duration 3). Therefore, to ensure that 
the Euclidean rhythm contains both binary and ternary cells, and no other durations, it follows 
that n must be between 2k and 3A;. 

Of course, n ot all a ksak rhythms are Euclidean. Consider the Bulgarian rhythm with i nterval 
sequence (3322) |Aro04l |. which is also the metric pattern of Indian Lady by Don Ellis [Kei9ll | . Here 
A; = 4 and n = 10, and £(4, 10) = [x • • x • x • • x • ] or (3232), a periodic rhythm. 
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The following Euclidean rhythms are authentic aksak: 

E(2,5) = [x • X • • ] = (23) (classical music, jazz, Greece, Macedonia, Namibia, Persia, Rwanda). 
^(3,7) = [x • X • X • • ] = (223) (Bulgaria, Greece, Sudan, Turkestan). 



^(4,11) = 


[x 


• X • • X • 






= (3332) (Southern India rhythm), (Serbian necklace). 


^(5,11) = 


[x 


X • X • X • 


X 


• •] 


= (22223) (classical music, Bulgaria, Northern India, Serbia) 


£^(5, 13) = 


[x 


• X • X • • 


X 


• X 


• ] = (32323) (Macedonia). 


S(6, 13) = 


[x 


X • X • X • 


X 


• X 


• ] = (222223) (Macedonia). 


E{7, 17) = 


[x 


• X • X • • 


X 


• X 


. X • X • ] = (3232322) (Macedonian necklace). 


Ei8, 17) = 


[x 


X • X • X • 


X 


• X 


X • X • • ] = (22222223) (Bulgaria). 


S(8, 19) = 


[x 


• X • X • X 




• X 


X • X • • X • ] = (32232232) (Bulgaria). 


S(9,23) = 


[x 


• X • X • • 


X 


• X 


•x-x - - x- x- •] = (323232323) (Bulgaria). 



The following Euclidean rhythms are quasi- aksak: 



S(4,9) = [x-x-x-x - •] = (2223) (Greece, Macedonia, Turkey, Zaire). 
S(7, 15) = [x • X • X • X • X • X • X • • ] = (2222223) (Bulgarian necklace). 

The following Euclidean rhythms are pseudo-aksak: 

E{3,8) = [x • • X • • X • ] = (332) (Central Africa, Greece, India, Latin America, West Africa, Su- 
dan). 

E{5, 12) = [x • • X • X • • X • X • ] = (32322) (Macedonia, South Africa). 

E{7, 16) = [x - - x- x- x- •x-x-x-] = (3223222) (Brazilian, Macedonian, West African neck- 
laces). 

E{7, 18) = [x • • X • X • • X • X • • X • X • • ] = (3232323) (Bulgaria). 
£;(9,22) = [x- - x- x- - x- x- - x- x- •x-x-] = (323232322) (Bulgarian necklace). 
£^(11,24) = [x- - x- x- x- x- x- .x-x-x-x-x-] = (32222322222) (Central African and Bul- 
garian necklaces). 

£(15,34) = [x- - x- x- x- x- - x- x- x- x- - x- x- x- x- •x-x-] = (322232223222322) (Bul- 
garian necklace). 



2.6 Drawing Digital Straight Lines 



Euclidean rhythms an d neck lace patterns also appear in the computer graphics literature on drawing 
digital straight lines KR04i |. The problem here consists of efficiently converting a mathematical 
straight line segment defined by the x and y integer coordinates of its endpoints, to an ordered 
sequence of pixels that most faithfully represents the given straight line segment. Figure [5] illustrates 
an example of a digital straight line (shaded pixels) determined by the two given endpoints p and 
q. All the pixels intersected by the segment (p, q) are shaded. If we follow either the lower or upper 
boundary of the shaded pixels from left to right we obtain the interval sequences (43333) or (33334), 
respectively. Note th at the upper pattern corresponds to £(5,16), a Bossa-N oya var iant. Indeed, 
Harris and Reingold HR04] show that the well-known Bresenham algorithm Bre65l | is described 
by the Euclidean algorithm. 



2.7 Calculating Leap Years in Calendar Design 

For thousands of years human beings have observed and measured the time it takes between two 
consecutive sunrises, and between tw o consecutive spring seasons. These measurements inspired 
different cultures to design calendars jAsc02l . IRD01[ | . Let Ty denote the duration of one revolution 
of the earth around the sun, more commonly known as a year. Let denote the duration of one 
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Figure 5: The shaded pixels form a digital straight line determined by the points p and q. 



complete rotation of the earth, more commonly known as a day. The values of Ty and are of 
course continually changing, because the universe is continually reconfiguring itself. However the 

ratio Ty/Td is approximately 365.242199 It is very convenient therefore to make a year last 365 

days. The problem that arises both for history and for predictions of the future, is that after a 

while the 0.242199 starts to contribute to a large error. One simple solution is to add one extra 

day every 4 years: the so-called Julian calendar. A day with one extra day is called a leap year. 

But this assumes that a year is 365.25 days long, which is still slightly greater than 365.242199 

So now we have an erro r in th e opposite direction albeit smaller. One solution to this problem is 
the Gregorian calendar |Sha94 . The Gregorian calendar defines a leap year as one divisible by 4, 



except not those divisible by 100, except not those divisible by 400. With this rule a year becomes 
365 + 1/4 - 1/100 + 1/400 = 365.2425 days long, not a bad approximation. 



Another solution is provided by the Jewish calendar which uses the idea of cycles Asc02l | . Here 
a regular year has 12 months and a leap year has 13 months. The cycle has 19 years including 7 
leap years. The 7 leap years must be distributed as evenly as possible in the cycle of 19. The cycle 
is assumed to start with Creation as year 1. If the year modulo 19 is one of 3, 6, 8, 11, 14, 17, or 
19, then it is a leap year. For example, the year 5765 = 303 -19-1-8 and so is a leap year. The year 
5766, which begins at sundown on the Gregorian date of October 3, 2005, is 5766 = 303x19 -|- 9, 
and is therefore not a leap year. Applying Bjorklund's algorithm to the integers 7 and 19 yields 
E{J, 19) = [x • • X • X • • X • • X • X • • X • • ]. If we start this rhythm at the 7th pulse we obtain 
the pattern [• - x- - x- x- - x- - x- - x- x], which describes precisely the leap year pattern 3, 
6, 8, 11, 14, 17, and 19 of the Jewish calendar. In this sense the Jewish calendar is an instance of 
a Euclidean necklace. 

2.8 Euclidean Strings 

In the study of t he combin atorics of words and sequences, there exists a family of strings called 
Euchdean strings jERSSOSl ]. In this section we explore the relationship between E uclidean s trings 



and Euclidean rhythms. We use the same terminology and notation introduced in [ERSS03 |. 

Let P = {po,pi, ...,pn-i) denote a string of non-negative integers. Let p{P) denote the right 
rotation of P by one position; that is, p{P) = {pn-i,Po,Pi, ■■■,Pn-2)- Let p'^{P) denote the right 
rotation of P by d positions. If P is considered as a cyclic string, a right rotation corresponds to a 
clockwise rotatio n. Fig ure [6| illustrates the p{P) operator with P equal to the Bembe bell-pattern 
of West Africa [T ou03] . Figure (Ha) shows the IB embe bell-pattern. Figure [6l[b) shows p{P), which 
is a hand-clapping pattern from W est Africa [Pre83i] . and Figure [6)[c) shows p'^{P), which is the 
Tambu rhythm of Curac ao Ros02[ | . 



Ellis et al. [ERSS03I ] define a string P = {pQ,pi, ...,pn~i) to be Euclidean if incrementing po 
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(a) 





Figure 6: Two right-rotations of tlie Bembe string: (a) the Bembe, (b) rotation by one unit, (c) rotation 
by seven units. 



by 1 and decrementing Pn-i by 1 yields a new string t{P) that is the rotation of P. In other 
words, P and t{P) are instances of the same necklace. Therefore, if we represent rhythms as 
binary sequences, Euclidean rhythms cannot be Euclidean strings because all Euclidean rhythms 
begin with a 'one'. Increasing pQ by one makes it a 'two', which is not a binary string. Therefore, 
to explore the relationship between Euclidean strings and Euclidean rhythms, we will represent 
rhythms by their clockwise distance sequences, which are also strings of nonnegative integers. As 
an example, consider ^(4,9) = [x-x-x-x - •] = (2223). Now r(2223) = (3222), which is a 
rotation of ^(4,9), and thus (2223) is a Euclidean string. Indeed, for P = £^(4,9), r(P) = p^{P). 
As a second example, consider the West African clapping-pattern shown in Figure [6)[b) given by 
P = (1221222). We have that t{P) = (2221221) = p^{P), the pattern shown in Figure[S^c), which 
also happens to be the mirror image of P about the (0, 6) axis. Therefore P is a Euclidean string. 
However, note that P is not a Euclidean rhythm. Nevertheless, P is a rotation of the Euclidean 
rhythm E{7, 12 ) = (2122 122). 

Ellis et al. [EE,SSn3t have many beautiful results about Euclidean strings. They show that 
Euclidean strings exist if, and only if, n and {po +pi + ...+Pn-i) are relatively prime numbers, and 
that when they exist they are unique. They also show how to construct Euclidean strings using an 
algorithm that has the same structure as the Euclidean algorithm. In addition they relate Euclide an 
strings to many other families of sequences studied in the combinatorics of words jASOilLotO^ . 

Let R{P) denote the reversal (or mirror image) of P; that is, R{P) = {pn^i,Pn-2, ■■■,Pi,Po)- 
Now we may determine which of the Euclidean rhythms used in world music listed in the preceding, 
are Euclidean strings or reverse Euclidean strings. The length of a Euclidean string is defined as the 
number of integers it has. This translates in the rhythm domain to the number of onsets a rhythm 
contains. Furthermore, strings of length one are Euclidean strings, trivially. Therefore all the trivial 
Euclidean rhythms with only one onset, such as £^(1,2) = [x • ] = (2), ^(1,3) = [x • • ] = (3), and 
S(l,4) = [x • • • ] = (4), etc., are both Euclidean strings as well as reverse Euclidean strings. In 
the lists that follow the Euclidean rhythms are shown in their box-notation format as well as in 
the clockwise distance sequence representation. The styles of music that use these rhythms is also 
included. Finally, if only a rotated version of the Euclidean rhythm is played, then it is still included 
in the list but referred to necklace. 

The following Euclidean rhythms are Euclidean strings: 

P(2,3) = [xx • ] = (12) (West Africa, Latin America, Nubia, Northern Canada). 

P(2, 5) = [x • X • • ] = (23) (classical music, jazz, Greece, Macedonia, Namibia, Persia, Rwanda), 

(authentic aksak). 
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E{3,4) = [xxx •] = (112) (Brazil, Bali rhythms), (Colombia, Greece, Spain, Persia, Trinidad 
necklaces) . 

£^(3,7) = [x • X • X • • ] = (223) (Bulgaria, Greece, Sudan, Turkestan), (authentic aksak). 
E{4:, 5) = [x X X X ■ ] = (1112) (Greece). 

i?(4, 9) = [x-x-x-x - •] = (2223) (Greece, Macedonia, Turkey, Zaire), (quasi-aksak). 

E{5, 6) = [x X X X X • ] = (11112) (Arab). 

E{5, 11) = [x • X • X • X • X • • ] = (22223) (classical music, Bulgaria, Northern India, Serbia), (au- 
thentic aksak). 

E(5, 16) = [x • • X • • X • • X • • X • • • • ] = (33334) (Brazihan, West African necklaces). 

E(6, 7) = [xxxxxx-] = (111112) (Greek necklace) 

£^(6, 13) = [x • X • X • X • X • X • • ] = (222223) (Macedonia), (authentic aksak). 
E(7,8) = [xxxxxxx •] = (1111112) (Libyan necklace). 

E(7, 15) = [x • X • X • X • X • X • X • • ] = (2222223) (Bulgarian necklace), (quasi-aksak). 
E(8, 17) = [x • X • X • X • X • X • X • X • • ] = (22222223) (Bulgaria), (authentic aksak). 

The following Euclidean rhythms are reverse Euclidean strings: 

E(3,5) = [x • X • x] = (221) (Korean, Rumanian, Persian necklaces). 

£^(3,8) = [x • • X • • X • ] = (332) (Central Africa, Greece, India, Latin America, West Africa, Su- 
dan), (pseudo- aksak). 

£(3, 11) = [x • • • X • • • X • • ] = (443) (North India). 
£(3, 14) = [x ■ ■ ■ • X • • • • X • • • ] = (554) (North India). 

£(4,7) = [x • X • X • x] = (2221) (Bulgaria). 

£(4, ll) = [x - - x- - x- •x-] = (3332) (Southern India rhythm), (Serbian necklace), (authentic 
aksak). 

£(4, 15) = [x ■ ■ ■ x • • • X • ■ ■ x ■ ■ ] = (4443) (North India). 

£;(5,7) = [x • XX • xx] = (21211) (Arab). 
£(5,9) = [x • X • X • X • x] = (22221) (Arab). 

£(5, 12) = [x • • X • X • • X • X • ] = (32322) (Macedonia, South Africa), (pseudo-aksak). 
E(7,9) = [x-xxx-xxx] = (2112111) (Greece). 
£(7,10) = [x-xx-xx-xx]= (2121211) (Turkey). 

E(7, 16) = [x - - x- x- x- •x-x-x-] = (3223222) (Brazihan, Macedonian, West African neck- 
laces), (pseudo-aksak). 

E(7, 17) = [x • • X ■ X • ■ X • X ■ ■ X ■ X ■ ] = (3232322) (Macedonian necklace), (authentic aksak). 
£(9,22) = [x- - x- x- - x- x- - x- x- .x-x-] = (323232322) (Bulgarian necklace), (pseudo- 
aksak). 

£(11, 12) = [x • X X X X X X X X X x] = (11111111112) (Oman necklace). 

£(11,24) = [x- - x- x- x- x- x- •x-x-x-x-x-] = (32222322222) (Central African and Bul- 
garian necklaces), (pseudo-aksak). 

The following Euclidean rhythms are neither Euclidean nor reverse Euclidean strings: 

£(5,8) = [x • XX • XX • ] = (21212) (Egypt, Korea, Latin America, West Africa). 

£(5, 13) = [x • • X • X • • X • X • • ] = (32323) (Macedonia), (authentic aksak). 

£(7,12) = [x-xx-x-xx-x-] = (2122122) (West Africa), (Central African, Nigerian, Sierra 

Leone necklaces). 

E(7, 18) = [x ■ ■ X ■ X ■ ■ X ■ X ■ ■ X ■ X ■ • ] = (3232323) (Bulgaria), (pseudo-aksak). 
E(8, 19) = [x • -x-x-x • - X - X - X - - X - ] = (32232232) (Bulgaria), (authentic aksak). 
£(9, 14) = [x-xx-xx-xx-xx-] = (212121212) (Algerian necklace). 



19 



E(9, 16) = [x-xx-x-x-xx-x-x-] = (212221222) (West and Central African, and Brazilian 
necklaces). 

£^(9,23) = [x- - x- x- - x- x- - x- x- - x- x- •] = (323232323) (Bulgaria), {authentic aksak). 
£(13,24) = [x-xx-x-x-x-x-xx-x-x-x-x-] = (2122222122222) (Central African neck- 
lace) . 

£(15,34) = [x- - x- x- x- x- - x- x- x- x- - x- x- x- x- •x-x-] = (322232223222322) 
(Bulgarian necklace), (pseudo-aksak). 

These three groups of Euclidean rhythms reveal a tantalizing pattern. The Euclidean rhythms 
that are favored in classical music and jazz are also Euclidean strings (the first group). Furthermore, 
this group is not popular in African music. The Euclidean rhythms that are neither Euclidean 
strings nor reverse Euclidean strings (group three) fall into two categories: those consisting of 
clockwise distances 1 and 2, and those consisting of clockwise distances 2 and 3. The latter group 
is used only in Bulgaria, and the former is used in Africa. Finally, the Euclidean rhythms that 
are reverse Euclidean strings (the second group) appear to have a much wider appeal. Finding 
musicological explanations for the preferences apparent in these mathematical properties raises 
interesting ethnomusicological quest ions. 

The Euclidean strings defined in determine another family of rhythms, many of which 



are also used in world music but are not necessarily Euclidean rhythms. For example, (1221222) 
is an Afro-Cuban bell pattern. Therefore it would be interesting to explore empirically the re- 
lation between Euclidean strings and world music rhythms, and to determine formally the exact 
mathematical relation between Euclidean rhythms and Euclidean strings. 



3 Definitions and Notation 

Before we begin the more technical part of the paper, we need to define some precise mathematical 
notation for describing rhythms. 

Let Z+ denote the set of positive integers. For k,n £ Z+, let gcd(A:,n) denote the greatest 
common divisor of k and n. If gcd(A:,n) = 1, we call k and n relatively prime. For integers a < b, 
let [a,b] = {a,a + l,a + 2, . . . , b}. 

Let C be a circle in the plane, and consider any two points x, y on C. The chordal distance 
between x and y, denoted by d{x,y), is the length of the line segment xy; that is, d{x,y) is the 
Euclidean distance between x and y. The clockwise distance from x to y, or of the ordered pair 
{x,y), is the length of the clockwise arc of C from x to y, and is denoted by d{x,y). Finally, the 
geodesic distance between x and y, denoted by d{x, y), is the length of the shortest arc of C between 
X and y; that is, d{x,y) = mm{d{x , y) , d{y , x)} . 

A rhythm of timespan n is a subset of {0, 1, . . . , n — 1}, representing the set of pulses that are 
onsets in each repetition. For clarity, we write the timespan n as a subscript after the subset: {. . . }„. 
Geometrically, if we locate n equally spaced points clockwise around a circle C„ of circumference n, 
then we can view a rhythm of timespan n as a subset of these n points. We consider an element of 
Cn to simultaneously be a point on the circle and an integer in {0, 1, . . . , n — 1}. 

The rotation of a rhythm R of timespan n by an integer A > is the rhythm {{i + A) mod n : 
i € R}n of the same timespan n. The scaling of a rhythm R of timespan n by an integer a > 1 is 
the rhythm {ai : i G R}an of timespan an. 

Let R = {ro, ri, . . . , rk-i}n be a rhythm of timespan n with k onsets sorted in clockwise order. 
Throughout this paper, an onset rj will mean (rj mod k) niod n. Observe that the clockwise distance 
d{ri,rj) = (rj — ri) mod n. This is the number of points on Cn that are contained in the clockwise 
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arc {ri,rj] and is also known as the chromatic length CD9l[. 

The geodesic distance multiset of a rhythm R is the multiset of all nonzero pairwise geodesic 
distances; that is, it is the multiset {'(i{ri,rj) : ri,rj ^ R,ri rj}. The geodesic distance multiset 
has cardinality (2) • The multiplicity of a distance d is the number of occurrences of d in the geodesic 
distance multiset. 

A rhythm is Erdos-deep if it has (exactly) one distance of multiplicity i, for each i G [l,k — 1]. 
Note that these multiplicities sum to X^t/ ^ ~ (2) ' "^hich is the cardinality of the geodesic distance 
multiset, and hence these distances are all the distances in the rhythm. Every geodesic distance is 
between and [n/2\. A rhythm is Winograd-deep if every two distances from {1,2, . . . , [§J} have 
different multiplicity. 

A shelling of an Erdos-deep rhythm R is an ordering si, §2; • • • > Sfc of the onsets in R such that 
R — {si, S2, ■ ■ ■ , Si} is an Erdos-deep rhythm for i = 0, 1, . . . , fc. (Every rhythm with at most two 
onsets is Erdos-deep.) 

The evenness of rhythm R is the sum of all inter-onset chordal distances in R; that is, 
d{ri,rj). 

0<i<j<k-l 

The clockwise distance sequence of R is the circular sequence (do, di, . . . , dfc_i) where di = 
d!{ri,ri^i) for all i £ [0, A; — 1]. Observe that each di G Z+ and Y^- di = n. 

Observation 1. There is a one-to-one relationship between rhythms with k onsets and timespan n 
and circular sequences (do, di, • • • , dk-i) where each di S Z,~^ and di = n. 



4 Even Rhythms 

In this section we first describe three algorithms that generate even rhythms. We then characterize 
rhythms with maximum evenness and show that, for given numbers of pulses and onsets, the three 
described algorithms generate the unique rhythm with maximum evenness. As mentioned in the 
introduction, the measure of evenness considered here is the pairwise sum of chordal dista nces. 



Th e even rhythms characterized in this section were studied by Clough and Meyerson CM85I . 
CM8fil ] for the case where the numbers of pulses an d onse ts are relatively prime. This was sub- 
sequently expanded upon by Clough and Douthett |CD9ll ]. We revisit these results and provide 
an additional connection to rhythm s (and scales) that are obtained from the Euclidean algorithm. 



Most of these results are stated in 
much more streamlined. 



CD9l|. However our proofs are new, and in many cases are 



4.1 Characterization 

We first present three algorithms for computing a rhythm with k onsets, timespan n, for any k < n, 
that possess large evenness. 



The first algorithm is by Clough and Douthett CD9l[: 



Algorithm CLOUGH-DoUTHETT(n, /;:) 
1. return {[f J : i G [0,A;- 1]} 



Because k <n, the rhythm output by CLOUGH-DouTHETT(n, A;) has k onsets as desired. 



The second algorithm is a geometric heuristic implicit in the work of Clough and Douthett CD9l[ | : 
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Algorithm SNAP(n, k) 

1. Let D be a set of k evenly spaced points on C„ such that D n C„ = 0. 

2. For each point x £ D, let x' be the first point in C„ clockwise from x. 

3. return {x' : x G D} 



Because k < n, the clockwise distance between consecutive points in D in the execution of 
SNAP(n, /c) is at least that of consecutive points in C„. Thus, x' ^ y' for distinct x,y £ D, so 
Snap returns a rhythm with k onsets as desired. 

The third algorithm is a recursive algorithm in the same mold as Euclid's algorithm for greatest 
common divisors. The algorithm uses the clockwise distance sequence notation described in the 
introduction. The resulting rhythm always defines the same necklace as the Euclidean rhythms 
from Section 12.31 that is, the only difference is a possible rotation. 



Algorithm EuCLlDEAN(n, A;) 

1. if A; evenly divides n then return (^, . . . , ^) 

k 



2. a <— n mod k 

3. {xi,X2, ■■■ ,Xa) 

4. return ([^J,. 



- Euclidean(/!;, a) 



X2—1 



As a simple example, consider k = 5 and n = 13. The sequence of calls to Euclidean (n, /c) 
follows the same pattern as the Euclid algorithm for greatest common divisors from Section [2.1^ 
except that it now stops one step earlier: (13, 5), (5, 3), (3, 2), (2, 1). At the base of the recursion, 
we have Euclidean(2, 1) = (2) = [x •]. At the next level up, we obtain Euclidean(3, 2) = 
(1, 2) = [x X • ]. Next we obtain Euclidean(5, 3) = (2; 1, 2) = [x • x x • ]. Finally, we obtain Eu- 
CLIDEAn(13, 5) = (2, 3; 3; 2, 3) = [x • X • • x • • x • x • • ]. (For comparison, the Euclidean rhythm 
from Section [2?2] is £'(5,13) = (2,3,2,3,3), a rotation by 5.) 

We now show that algorithm EuCLiDEAN(n, /c) outputs a circular sequence of k integers that 
sum to n (which is thus the clockwise distance sequence of a rhythm with k onsets and timespan n). 
We proceed by induction on k. If k evenly divides n, then the claim clearly holds. Otherwise a 



(= n mod k) > 0, and by induction Yli= 
and sums to 



k. Thus the sequence that is output has k terms 







n 1 




+ 




k 


.k\ 



i=l 



n 



1) = a - +{k-a) 



a(l+ ^ ^ + {k-a) 
a + k 



n 



The following theorem is one of the main contributions of this paper. 

Theorem 4.1. Let n > k > 2 be integers. The following are equivalent for a rhythm R = 
{ro,ri, . . . ,rk-i}n with k onsets and timespan n: 

(A) R has maximum evenness (sum of pairwise inter-onset chordal distances), 
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(B) R is a rotation of the CLOUGH-DoUTHETT(n, k) rhythm, 

(C) R is a rotation of the SNAP(n, /c) rhythm, 

(D) R is a rotation of the Euclidean (n, /c) rhythm, 

(*) for all £ S [1, k] and i £ [0,k — 1], the ordered pair (r^, rj+^) has clockwise distance (t{ri, r^j^i) € 

Moreover, up to a rotation, there is a unique rhythm that satisfies these conditions. 

Note that the evenness of a rhythm equals the evenness of the same rhythm played backwards. 
Thus, if R is the unique rhythm with maximum evenness, then R is the same rhythm as R played 
backwards (up to a rotation). 

The proof of Theorem 14.11 proceeds as follows. In Section 14.21 we prove that each of the three 
algorithms produces a rhythm that satisfies property (*). Then in Section [4]3] we prove that there is 
a unique rhythm that satisfies property (*). Thus the three algorithms produce the same rhythm, 
up to rotation. Finally in Section 14.41 we prove that the unique rhythm that satisfies property (*) 
maximizes evenness. 



4.2 Properties of the Algorithms 



We now prove that each of the algorithms has property (★). Clough and Douthett CD9ll | proved 
the following. 

Proof (B) =^ (★). Say i? = {ro, n, . . . , rfc_i}„ is the CLOUGH-DoUTHETT(n, A;) rhythm. Consider 
an ordered pair (ri,ri^i) of onsets in R. Let pi = in mod k and let p£ = in mod k. By symmetry 
we can suppose that rj ^ r^j.j.^^ mod k- 

Then the clockwise distance d{ri,ri^e) is 
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which is L^J or [^], because L^i+Rj g {o,l}. 



□ 



A similar proof shows that the rhythm {\^~\ : « G [0, — 1]} satisfies property (★). Observe 
that (*) is equivalent to the following property. 



(:*n*r) if {do, di, . . . , dk-i) Is the clockwise distance sequence of R, then for all £ G [l, ^1, the sum of 
any £ consecutive elements in (do, di, . . . , dfc-i) equals [^] or [^J . 

Proof (C) (^)- Let {d^, di, . . . , d^-i) be the clockwise distance sequence of the rhythm deter- 
mined by SNAp(n, k). For the sake of contradiction, suppose that for some i G [1, k], the sum of i 
consecutive elements in {dQ,di, . . . ,dk-i) is greater than [^]. The case in which the sum is less 
than [^J is analogous. We can assume that these £ consecutive elements are {do,di, . . . ,d£^i). 
Using the notation defined in the statement of the algorithm, let xq, xi, . . . ,Xi be the points in D 
such that d(x-,x-^_i) = di for all i € [0,£- 1]. Thus d{x[, x'^^-^^) > \^] + 1. Now d{xe+l,x'^_^_-^) < 1. 
Thus 'ci{x[, X£^i) > > which implies that cl{xi,X£^i) > This contradicts the fact that 
the points in D were evenly spaced around C„ in the first step of the algorithm. 

□ 
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Proof (D) ^ We proceed by induction on /c. Let i? = Euclidean (n, A;). If A; evenly divides n, 

then R = which satisfies (D). Otherwise, let a = n mod k and let {xi,X2, ■ ■ ■ , Xa) = 

Euclidean(A;, a). By induction, for all i £ the sum of any i consecutive elements in 

{xi,X2, ■ ■ ■ ,Xa) equals [^J or [^]. Let 5 be a sequence of m consecutive elements in R. By 
construction, for some 1 < i < j < a, and for some < s < Xj — 1 and < t < x.- — 1, we have 



mn 
k 



It remains to prove that |_^J < 5* < p 

We first prove that J2 ^ — |_^J • assume the worst case for ^ 5 to be minimal, which 

is when s = Xj — 1 and t = Xj — 1. Thus by induction. 
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Now we prove that ^ S* < |_^J • We can assume the worst case for ^ to be maximal, which 
is when s = and t = 0. Thus by induction. 



i-i 

m — 1 = Xa > 

a=i+l 



l)k 



Hence 



am a 

~T-k 



(j 



DA; 



a a f ii — i — ^)k — a + \\ a 1 



Thus [^] >j-i and 



i < m 


n 










n 


am^ 








n 






'mn' 






A; 




m 




+ - 












k 




-k. 






-k- 















□ 



4.3 Uniqueness 

In this section we prove that there is a unique rhythm satisfying the conditions in Theorem 14. 1[ The 
following well-known number-theoretic lemmas will be useful. Two integers x and y are inverses 
modulo m if xy = 1 (mod m). 

Lemma 4.2 ( [StiOd . page 55]). An integer x has an inverse modulo m if and only if x and m are 
relatively prime. Moreover, if x has an inverse modulo m, then it has an inverse y £ [l,m — 1]. 
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Lemma 4.3. Ifx and m are relatively prime, then ix ^ jx (mod m) for all distinct i, j G [0, m — 1] . 

Proof. Suppose that ix = jx (mod m) for some i,j £ [0, m — 1]. By Lemma \i.2\ x has an inverse 
modulo m. Thus i = j (mod m), and i = j because i,j £ [0,m — 1]. □ 

Lemma 4.4. For all relatively prime integers n and k with 2 < k < n, there is an integer i E 
[1, A; — 1] such that: 

(a) in = 1 (mod k), 

(b) a ^ ji (mod k) for all distinct i,j £ [0,k — I], and 

(c) «LtJ ^ jlT\ n) for all distinct i,j G [0, - 1]. 

Proof. By Lemma [4.21 with x = n and m = k, n has an inverse ^ modulo A;. This proves (a). Thus k 
and £ are relative prime by Lemma 14.21 with x = i and m = k. Hence (b) follows from Lemma 14.31 
Let t = [^J. Then £n = kt + 1. By Lemma 14.31 with m = n and x = t (and because k < n), to 
prove (c) it suffices to show that t and n are relatively prime. Let g = gcd(t,n). Thus = A;l + 1 
Because ^ and | are integers, ^ is an integer and g = 1. This proves (c). □ 

The following theorem is the main result of this section. 

Theorem 4.5. For all integers n and k with 2 < k < n, there is a unique rhythm with k onsets 
and timespan n that satisfies property (★), up to a rotation. 

Proof. Let R = {ro,ri, . . . ,r'fc„i}„ be a fc-onset rhythm that satisfies (*). Recall that the index of 
an onset is taken modulo k, and that the value of an onset is taken modulo n. That is, ri = x 
means that rj mod k = x mod n. 

Let g = gcd(?7,, k). We consider three cases for the value of g. 

Case 1. g = k: Because R satisfies property (*) for £ = 1, every ordered pair (ri,rj+i) has 
clockwise distance By a rotation of R we can assume that ro = 0. Thus fi = ^ for all 
i £ [0,k — 1]. Hence R is uniquely determined in this case. 

Case 2. g = 1 (see Figure [7]): By Lemma 14.4( a). there is an integer £ £ [1,A; — 1] such that 
£n = 1 (mod k). Thus in = {k — l)LxJ + Txl • Hence, of the k ordered pairs (rj,rj_|_^) of onsets, 
k — 1 have clockwise distance [^J and one has clockwise distance \^~\ . By a rotation of R we 
can assume that tq = and = n — \^~\. Thus r^ = i[^\ for all i £ [0,k — 1]; that is, 

'>^{ii) mod k = (^LxJ) ™od n. By Lemma [4.4( b) and (c), this defines the k distinct onsets of R. 
Hence R is uniquely determined in this case. 

Case 3. g £ [2,k - 1] (see Figure E]): Let k' = ^ and let n' = |. Observe that both k' and 

n' are integers. Because R satisfies (*) and = L^J = ^'i have 'd[{ri,rij^ki) = n! for all 

i G [0,/c - 1]. Thus 

rifc/+j = in + rj (1) 

for ah i£ [0, g - 1] and j £ [0, n' - 1]. 

Now gcd(n'. A;') = 1 by the maximality of g. By Lemma [4.4( a). there is an integer i! £ [1, k' — 1] 
such that fn' = 1 (mod k'). Thus f n' = (A;'- 1) [^J + [^] , implying fn = {k-g)\_^\^g\^^'\. 
Hence, of the k ordered pairs (rj,rj_|_£/) of onsets, k — g have clockwise distance L^F"J and g have 
clockwise distance [^p^] • By a rotation of R we can assume that tq = and r^/ = [^p-] • By 
Equation ([1]) with j = and j = we have 



Tjfc/ = iri and rik'+v = + 



(2) 
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Figure 7: Here we illustrate Case 2 with n = 12 and k — 7. Thus £ = 3 because 3 x 12 = 1 (mod 7). We 
have \^~\ = 6 and [^J =5. By a rotation we can assume that tq — and Vk-i = ^4 = 6 (the darker dots). 
Then as shown by the arrows, the positions of the other onsets are implied. 



for all i G [0,g — !]• This accounts for the g ordered pairs (rj,rj_|_^/) with clockwise distance [^p-]- 
The other k — g ordered pairs {ri,ri^£i) have clockwise distance Lt?^J- Define 



-Lq = and L 



+ (j - 1) 



for all i £ [1, k' -1] . 



Thus by Equation ^ 



nk'+jP = in + Lj 

for all i G [0,g — 1] and j G [0, k' — 1]; that is, r(^if^/_^_j£/^ ^^^^ ^ = {in' + Lj) mod n. 






11 t 



Figure 8: Here we illustrate Case 3 with n = 15 and k ^ 9. Thus g = 3, n' = 5 and k' = 3. We have 
e' = 2 because 2x5=1 (mod 3). Thus = 4 and [^J = 3. We have = 0, Li = 4 and L2 = 7. A 

rotation fixes the first .g = 3 onsets (the darker or blue dots). As shown by the arrows, these onsets imply 
the positions of the next three onsets (medium or green dots), which in turn imply the positions of the final 
three onsets (the light or yellow dots). 

To conclude that R is uniquely determined, we must show that over the range i £ [0, g — 1] and 
j G [0, k' — 1], the numbers ik' + ji' are distinct modulo k, and the numbers in' + Lj are distinct 
modulo n. 
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First we show that the numbers ik' + ji' are distinct modulo k. Suppose that 

ik' + jl' = pk' + ji' (mod A;) 



(3) 



for some i,p € [0,g — 1] and j, q G [0, k' — 1]. Because k = k' ■ g, we can write {ik' + ji') mod k as 
a multiple of k' plus a residue modulo k'. In particular, 

(ifc' mod k = k' ({i+ [^J) mod 9) + {ji' mod A;') . 



Thus Equation ^ implies that 

^' ((^ + Lp^J ) mod 5) + (i/ mod A:') = k' (^{p + [^J ) mod 5) + (g/ mod k') . (4) 

Hence ji' = qi' (mod A;'). Thus j = q hy Lemma l4.4r c). By substituting j = q into Equation (jj]), 
it follows that i = p (mod g). Thus i = p because G [0,g — 1]. This proves that the numbers 
ik' + ji' are distinct modulo k. 

Now we show that the numbers in' + Lj are distinct modulo n. The proof is similar to the 
above proof that the numbers ik' + ji' are distinct modulo k. 

Suppose that 

in' + Lj = pn + Lg (mod n) (5) 

for some i,p £ [0,g — 1] and j, q £ [0, A;' — 1]. Because n = n' ■ g, we can write {in' + Lj) mod n as 
a multiple of n' plus a residue modulo n' . In particular, 



(in' + Lj) mod n = n' (^(i + ^ ) mod gj + {Lj mod n') 
Thus Equation ([5]) implies that 



n 



mod 5 ) + {Lj mod n') = n' ( (p + 



mod g I + (Lg mod n) . 



(6) 



Hence Lj = Lq (mod n'). We claim that j = q. If j = then Lj = 0, implying Lq = and q = 0. 
Now assume that j,q > 1. In this case, Lj = iL^p"J + 1 and Lq = q[^-^\ + 1- Thus 



i'n' 



i'n' 



(mod A;') . 

Hence j = q hy Lemma [4.4f c). By substituting j = q into Equation ([6]), it follows that i = p 
(mod g). Thus i = p because i,p G [0,g — 1]. This proves that the numbers in' + Lj are distinct 
modulo n. 

Therefore R is uniquely determined. □ 

We have shown that each of the three algorithms generates a rhythm with property (★), and 
that there is a unique rhythm with property (*). Thus all of the algorithms produce the same 
rhythm, up to rotation. It remains to prove that this rhythm has maximum evenness. 

4.4 Rhythms with Maximum Evenness 

We start with a technical lemma. Let v, w be points at geodesic distance d on a circle C. Obviously 
d{v, w) is a function of d, independent of v and w. Let /(C, d) = d{v, w). 
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Lemma 4.6. For all geodesic lengths x < d on a circle C , we have f{C,x)+f(C,d—x) <2-/(C, 
with equality only if d = 2x. 

Proof. We can assume that C is a unit circle. Consider the isosceles triangle formed by the center 
of C and a geodesic of length d (< vr). We have ^/(C, d) = sin |. Thus /(C, d) = 2 sin ^. Thus our 
claim is equivalent to sinx + sin((i — x) < 2sin | for all x < d (< 7r/2). In the range < x < d, sinx 
is increasing, and sm{d — x) is decreasing at the opposite rate. Thus sina; + sin(d — x) is maximized 
when X = d — X. That is, when d = 2x. The result follows. □ 

For a rhythm R = {j'oi ^ii • • • > '"fc-ilni for each £ £ [l,k], let S{R,i) be the sum of chordal 
distances taken over all ordered pairs (r"j,rj_|_^) in R. That is, let S{R,i) = "^iZo 
Property (A) says that R maximizes X]^=i S{R,i). Before we characterize rhythms that maximize 
the sum of S{R,i), we first concentrate on rhythms that maximize S{R,i) for each particular 
value of i. Let D{R,i) be the multiset of clockwise distances {d{ri,ri^i) : i G [0,k — 1]}. Then 
S{RJ) is determined by D{R,£). In particular, S{R,e) = E{/(Cn,d) : d G D{R,l)} (where 
{f{Cn,d) : d G D{R,e)} is a multiset). 

Lemma 4.7. Let 1 < £ < k < n be integers. A k-onset rhythm R = {ro,ri, . . . ,rfc_i}„ maximizes 
S{R,i) if and only if \^{ri,ri^£) — ^{rj,rj^i)\ < 1 for all i,j £ [0,k — 1]. 

Proof. Suppose that R = {tq, n, . . . , rfe_i}„ maximizes S{R,i). Let di = d{ri,ri^£) for all i G 
[0, k — 1]. Suppose on the contrary that dp > dq + 2 for some p,q £ [0,k — 1]. We can assume 
that q < p, dp = dg + 2, and di = dq + 1 for all i G [q + l,p — 1]. Define = rj + 1 for all 
i G [q + and define r'- = ri for all other i. Let R' be the rhythm {r'Q,r^, . . . ,r^_]^}n- Thus 

D{RJ)\D{R'J) = {dp,dq} and D{R',e)\D{R,£) = {dp-l,dq+l}. Nowdp-1 = dq+1 = l{dp+dq). 
By Lemma Ha f{Cn,dp) + f{Cn,dq) < 2 ■ /(C„,i(dp + dq). Thus S{R,£) < S{R',£), which 
contradicts the maximality of S{R,i). 

For the converse, let Rhe a rhythm such that \(^{ri,ri^£) — Si^r^, rj_|_£)| < 1 for all i, j G [0, A; — 1]. 
Suppose on the contrary that R does not maximize S{R, I). Thus some rhythm T = {to, ti, . . . , tfc-i) 
maximizes S{T,£) and T ^ R. Hence D{T,i) / D{R,i). Because Y.D{R,£) = Y.D{T,£) (= £n), 
we have (^{ti,ti^£) — (t{tj,tj^i) > 2 for some i,j G [0, — 1]. As we have already proved, this implies 
that T does not maximize S{T,£). This contradiction proves that R maximizes S{R,£). □ 

Because ^i=Q d{ri,ri^i) = £n for any rhythm with k onsets and timespan n. Lemma 14.71 can 
be restated as follows. 

Corollary 4.8. Let l<£<k<nbe integers. A k-onset rhythm R = {ro,ri, . . . , maximizes 
S{R,i) if and only ifd{r,,r,+,) G {[f 1, L^J} V all i £ [0,fc- 1]. □ 

Proof (*) (A). If (★) holds for some rhythm i?, then by Corollarv 14.81 R maximizes S{R,£) for 
every £. Thus R maximizes S{R,£). □ 

Proof (A) ^ (★) . By Theorem 14. 5^ there is a unique rhythm R that satisfies property (★) . Let R 
denote the unique rhythm that satisfies property (*). Suppose on the contrary that there is a rhythm 
T = (to;^i) • • • lik-i) with property (A) but R ^ T. Thus there exists an ordered pair in 
T with clockwise distance d{ti,ti+i) {[f J, fxl}- Corollary HSl S{T,£) < S{R,£). Because 
T has property (A), ELi'5(r,£) > Y^i=iS{R,£). Thus for some £' we have S{T,i') > S{R,e'). 
But this is a contradiction, because S{R,£') > S{T,£') by Corollarv 14.81 □ 
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This completes the proof of Theorem I4.1[ Note that Theorem 14.11 holds for metrics other than 
the pairwise sum of all chordal distances. The only property that is needed is Lemma 14.61 For 
example, the metric "pairwise sum of the squares of all geodesic distances" satisfies Lemma 14.61 
and thus Theorem 14.11 holds in this setting. 



5 Deep Rhythms 



Winograd [Win66| | . and independently Clough et al. CEK99| . characterize all Winograd-deep scales: 



up to rotation, they are the scales that can be generated by the first [ri/2j or [n/2j +1 multiples 
(modulo n) of a value that is relatively prime to n, plus one exceptional scale {0, l,2,4}e. In 
this section, we prove a similar (but more general) characterization of Erdos-deep rhythms: up to 
rotation and scaling, they are the rhythms generable as the first k multiples (modulo n) of a value 
that is relatively prime to n, plus the same exceptional rhythm {0, 1,2,4}6. The key difference is 
that the number of onsets k is now a free parameter, instead of being forced to be either [n/2\ or 
[n/2j + 1. Our proof follows Winograd's, but differs in one case (the second case of Theorem 15. 3p . 

We later prove that every Erdos-deep rhythm has a shelling and that maximally even rhythms 
with n and k relatively prime are Erdos-deep. 



5.1 Characterization of Deep Rhythms 

Our characterization of Erdos-deep rhythms is in terms of two families of rhythms. The main 
rhythm family consists of the generated rhythms Dk^n,m = {^"^ mod n : i = 0,1, . . . ,k — 1}„ of 
timespan n, for certain values of k, n, and m. The one exceptional rhythm is -F = {0, l,2,4}e of 
timespan 6. 

Fact 5.1. F is Erdos-deep. 

Lemma 5.2. If k < [ra/2j -|- 1 and m and n are relatively prime, then Di^^n,m is Erdos-deep. 

Proof. The multiset of clockwise distances in Z?fc,n,m is {{jm — im) mod n : i < = {{j — i)m mod 
n : i < j}. There are k—p choices of i and j such that j — i = p, so there are exactly p occurrences of 
the clockwise distance (pm) mod n in the multiset. Each of these clockwise distances corresponds 
to a geodesic distance — either (pm) mod n or {—pm) mod n, whichever is smaller (at most n/2). 
We claim that these geodesic distances are all distinct. Then the multiplicity of each geodesic 
distance {±pm) mod n is exactly p, establishing that the rhythm is Erdos-deep. 

For two geodesic distances to be equal, we must have zizpm = zizqm (mod n) for some (possibly 
different) choices for the it symbols, and for some p ^ q. By (possibly) multiplying both sides 
by —1, we obtain two cases: (1) pm = qm (mod n) and (2) pm = —qm (mod n). Because m is 
relatively prime to n, by Lemma 14.21 "t- has a multiplicative inverse modulo n. Multiplying both 
sides of the congruence by this inverse, we obtain {1) p = q (mod n) and (2) p = —q (mod n). 
Because 0<'j<j<fc< [n/2\ -\- 1, we have 0<p = j — i< [n/2\ + 1, and similarly for q: 
< p,q < \n/2\. Thus, the first case of p = q (mod n) can happen only when p = q, and the 
second case of p + g = (mod n) can happen only when p = q = or when p = q = n/2. Either 
case contradicts that p ^ q. Therefore the geodesic distances arising from different values of p are 
indeed distinct, proving the lemma. □ 

We now state and prove our characterization of Erdos-deep rhythms, which is up to rotation 
and scaling. Rotation preserves the geodesic distance multiset and therefore Erdos-deepness (and 



29 



Winograd-deepness). Scaling maps each geodesic distance d to ad, and thus preserves multiphcities 
and therefore Erdos-deepness (but not Winograd-deepness). 

Theorem 5.3. A rhythm is Erdos-deep if and only if it is a rotation of a scaling of either the 
rhythm F or the rhythm Dk^n,m for some k,n,m with k < [n/2\ +1, 1 <m < \n/2\, and m and 
n are relatively prime. 

Proof. Because a rotation of a scaling of an Erdos-deep rhythm is Erdos-deep, the "if" direction of 
the theorem follows from Fact 15.11 and Lemma 15.21 

Consider an Erdos-deep rhythm R with k onsets. By the definition of Erdos-deepness, R 
has one nonzero geodesic distance with multiplicity i for each i = 1, 2, . . . , /c — 1. Let m be the 
geodesic distance with multiplicity k — 1. Because m is a geodesic distance, 1 < m < [n/2j. Also, 
k < \n/2\ + 1 (for any Erdos-deep rhythm R), because all nonzero geodesic distances are between 
1 and [n/2j and therefore at most [n/2j nonzero geodesic distances occur. Thus k and m are 
suitable parameter choices for Dk^n,m- 

Consider the graph Gm = Em) with vertices corresponding to onsets in R and with an edge 
between two onsets of geodesic distance m. By the definition of geodesic distance, every vertex i in 
Gm has degree at most 2: the only onsets at geodesic distance exactly m from i are {i — m) mod n 
and (i + m) mod n. Thus, the graph Gm is a disjoint union of paths and cycles. The number of 
edges in Gm is the multiplicity of m, which we supposed was k — 1, which is 1 less than the number 
of vertices in Gm- Thus, the graph Gm consists of exactly one path and any number of cycles. 

The cycles of Gm have a special structure because they correspond to subgroups generated 
by single elements in the cyclic group (Z/(n),-|-). Namely, the onsets corresponding to vertices 
of a cycle in Gm form a regular (n/a)-gon, with a geodesic distance of a = gcd(m, n) between 
consecutive onsets, (a is called the index of the subgroup generated by m.) In particular, every 
cycle in Gm has the same length r = n/a. Because Gm is a simple graph, every cycle must have at 
least 3 vertices, so r > 3. 

The proof partitions into four cases depending on the length of the path and on how many 
cycles the graph Gm has. The first two cases will turn out to be impossible; the third case will lead 
to a rotation of a scaling of rhythm F; and the fourth case will lead to a rotation of a scaling of 
rhythm Dk,n,m- 

First suppose that the graph Gm consists of a path of length at least 1 and at least one cycle. 
We show that this case is impossible because the rhythm R can have no geodesic distance with 
multiplicity 1. Suppose that there is a geodesic distance with multiplicity 1, say between onsets ii 
and 22- If i is a vertex of a cycle, then both (i + m) mod n and (i — m) mod n are onsets in R. If 
i is a vertex of the path, then one or two of these are onsets in R, with the case of one occurring 
only at the endpoints of the path. If (ii -|- m) mod n and (z2 + m) mod n were both onsets in R, or 
(ii — m) mod n and {i2 — m) mod n were both onsets in R, then we would have another occurrence of 
the geodesic distance between ii and i2, contradicting that this geodesic distance has multiplicity 1. 
Thus, ii and 12 must be opposite endpoints of the path. If the path has length i, then the clockwise 
distance between ii and Z2 is (Im) mod n. This clockwise distance (and hence the corresponding 
geodesic distance) appears in every cycle, of which there is at least one, so the geodesic distance 
has multiplicity more than 1, a contradiction. Therefore this case is impossible. 

Second suppose that the graph Gm consists of a path of length and at least two cycles. We 
show that this case is impossible because the rhythm R has two geodesic distances with the same 
multiplicity. Pick any two cycles G and C", and let d be the smallest positive clockwise distance 
from a vertex of C to a vertex of C'. Thus z is a vertex of C if and only if (i -|- d) mod n is a vertex 
of G'. Because the cycles are disjoint, d < a. Because r > 3, d < n/3, so clockwise distances of 
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d are also geodesic distances of d. The number of occurrences of geodesic distance d between a 
vertex of C and a vertex of C is either r or 2r, the case of 2r arising when d = a/2 (that is, C is a 
"half-rotation" of C). The number of occurrences of geodesic distance d' = min{d-|-m, n— (d + m)} 
is the same — either r or 2r, in the same cases. (Note that d < a < n — m, so d + m < n, so 
the definition of d' correctly captures a geodesic distance modulo n.) The same is true of geodesic 
distance d" = mm{d — m,n — {d — m)}. If other pairs of cycles have the same smallest positive 
clockwise distance d, then the number of occurrences of d, d' , and d" between those cycles are also 
equal. Because the cycles are disjoint, geodesic distance d and thus d + m and d — m cannot be 
{pm) mod n for any p, so these geodesic distances cannot occur between two vertices of the same 
cycle. Finally, the sole vertex x of the path has geodesic distance d to onset i (which must be a 
vertex of some cycle) if and only if x has geodesic distance d' to onset {i + m) mod n (which must 
be a vertex of the same cycle) if and only if x has geodesic distance d" to onset {i — m) mod n 
(which also must be a vertex of the same cycle). Therefore the multiplicities of geodesic distances 
d, d', and d" must be equal. Because R is Erdos-deep, we must have d = d' = d" . To have d = d', 
either d = d-\-raoid = n — (d + m), but the first case is impossible because d > by nonoverlap of 
cycles, so 2d + m = n. Similarly, to have d = d" , we must have 2d — m = n. Subtracting these two 
equations, we obtain that 2?7i = 0, contradicting that m > 0. Therefore this case is also impossible. 

Third suppose that the graph Gm consists of a path of length and exactly one cycle. We show 
that this case forces -R to be a rotation of a scaling of rhythm F because otherwise two geodesic 
distances m and m' have the same multiplicity. The number of occurrences of geodesic distance m 
in the cycle is precisely the length r of the cycle. Similarly, the number of occurrences of geodesic 
distance m' = min{2m, n — 2m} in the cycle is r. The sole vertex x on the path cannot have 
geodesic distance m or m' to any other onset (a vertex of the cycle) because then x would then be 
on the cycle. Therefore the multiplicities of geodesic distances m and m' must be equal. Because 
R is Erdos-deep, m must equal m', which implies that either m = 2m or m = n — 2m. The first 
case is impossible because m > 0. In the second case, 3m = n, that is, m = |n. Therefore, the 
cycle has r = 3 vertices, say at A, A -|- |n, A -|- |n. The fourth and final onset x must be midway 
between two of these three onsets, because otherwise its geodesic distance to the three vertices are 
all distinct and therefore unique. No matter where x is so placed, the rhythm R is a rotation by 
A -|- c^n (for some c G {0, 1, 2}) of a scaling by n/6 of the rhythm F. 

Finally suppose that Gm has no cycles, and consists solely of a path. We show that this case 
forces to be a rotation of a scaling of a rhythm Dk,n',m' with 1 <m' < [n'/2j and with m' and n' 
relatively prime. Let h be the onset such that (b — m) mod n is not an onset (the "beginning" vertex 
of the path). Consider rotating R by —i so that is an onset in the resulting rhythm R — i. The 
vertices of the path in i? — i form a subset of the subgroup of the cyclic group (Z/ (n), -|-) generated 
by the element m. Therefore the rhythm R — i = D}^^n,m = {{im) mod n : i = 0, l,...,fc — 1}„ is 
a scaling by a of the rhythm Di^^n/a,m/a = {{im/a) mod (n/a) : i = 0,1, . . . , k — 1}„. The rhythm 

.n/a,m/a ^as an appropriate value for the third argument: m/a and n/a are relatively prime 
(a = gcd(m,n)) and 1 < m/a < \n/2\/a < [{n/a)/2\. Also, k < [{n/a)/2\ + 1 because the only 
occurring geodesic distances are multiples of a and therefore the number A: — 1 of distinct geodesic 
distances is at most [{n/a)/2\. Therefore R is a rotation by i of a scaling by a of Dj^ n/a,m/a with 
appropriate values of the arguments. □ 

Corollary 5.4. A rhythm is Erdos-deep if and only if it is a rotation of a scaling of the rhythm F or 
it is a rotation of a rhythm Di^ n,m for some k, n, m satisfying k < \n/2g\ + 1 where g = gcd(m, n). 

Proof. First we show that any Erdos-deep rhythm has one of the two forms in the corollary. By 
Theorem 15. 3| there are two fiavors of Erdos-deep rhythms, and the corollary directly handles 
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rotations of scalings of F. Thus it suffices to consider a rhythm R that is a rotation by A of a scahng 
by a of Dk,n,m where k < [n/2j +1, 1 < m < [n/2\, and m and n are relatively prime. Equivalently, 
i? is a rotation by A of Dk,n',m' where n' = an and m' = am. Now g = gcd(n', m!) = a, so n' / g = n. 
Hence, k < \n' /2g\ + 1 as desired. Thus we have rewritten R in the desired form. 

It remains to show that every rhythm in one of the two forms in the corollary is Erdos-deep. 
Again, rotations of scalings of F are handled directly by Theorem 15. 31 So consider a rotation 
of Dk^nm where k < \n/2g\ + 1. The value of m matters only modulo n, so we assume that 
< m < n — 1. 

First we show that, if [n/2j +1 < m < n — 1, then D}^ n,m can be rewritten as a rotation of 
the rhythm Dk,n,m' where m! = n — m < \n/2\. By reversing the order in which we list the onsets 
in Dk,n,m = mod n : i = 0, 1, . . . , /c — we can write Dk,n,m = {{k — 1 — i)m mod n : i = 
0,1, ... ,k— l}n. Now consider rotating the rhythm Dkn,n-m = {i{n — m) mod n : i = . . . ,k — 
l}n by {k — l)m. We obtain the rhythm {[i (n — m) + {k — 1) m] mod n : i = 0, 1, . . . , — 1}„ = 
{[{k — l—i) m+in] mod n : i = 0, 1, . . . , k—\}n = {(^—1 — ^) "t- mod n : i = 0, 1, . . . , fc— l}n = F)k^n,m 
as desired. 

Thus it suffices to consider rotations of Dk^n,m where 1 < m < [n/2\ and k < [n/2g\. The 
rhythm Dk,n',m'^ where n' = n/g and m' = m/g, is Erdos-deep by Theorem 15.31 because n' and m' 
are relatively prime, k < \n' /2\ + 1, and 1 <m' < [n'/2j . But D}^^n,m is the scaling of Dk^n',m' by 
the integer g, so Dk^n,m is also Erdos-deep. □ 

An interesting consequence of this characterization is the following: 

Corollary 5.5. Every Erdos-deep rhythm has a shelling. 

Proof. If the Erdos-deep rhythm is Dk^n,m, we can remove the last onset from the path, resulting 
in Dk-i^n,m, and repeat until we obtain the empty rhythm -Do,n,m- At all times, k remains at most 
[n/2\ +1 (assuming it was originally) and m remains between 1 and [n/2j and relatively prime to n. 
On the other hand, F = {0, 1, 2, 4}6 has the shelling 4, 2, 1, because {0, 1, 2}e is Erdos-deep. □ 



5.2 Connection Between Deep and Even Rhythms 

A connectio n between maximally even scales and Winograd-deep scales is shown by Clough et 



al. [CEK99|. They define a diatonic scale to be a maximally even scale with k = {n -\- 2)/2 and n 
a multiple of 4. They show that diatonic scales are Winograd-deep. We now prove a similar result 
for Erdos-deep rhythms. 

Lemma 5.6. A rhythm R of maximum evenness satsifying k < [n/2\ +1 is Erdos-deep if and only 
if k and n are relatively prime. 

Proof. Recall that by property (*) one of the unique characterizations of an even rhythm of max- 
imum evenness can be stated as follows. For all 1 < £ < A;, and for every ordered pair (rj,rj_|_£) of 
onsets in R, the clockwise distance 'S{ri,rij^i) G {LxJ' Txli' 

For the case in which k and n are relatively prime, by Lemma 14.2^ there exists a value i < k 
such that in = 1 (mod k). Thus we can write in = k[in/k\ -\- 1. Let m = [in/k\. Now consider 
the set {im mod n : i = 0,1, . . . , k — 1}„. By Lemma H^ c) . we get k distinct values, so R can be 
realized as Dk^n,m = {i^ mod n : i = 0,1, . . . ,k — 1}„. Thus, by Lemma [5?2l R is Erdos-deep. 

Observe that F = {0, l,2,4}e does not maximize evenness because 3'(0, 2) = 2 and d[{2, 0) = 4 
yet i = 2. Hence, any rhythm that maximizes evenness and that is deep must also be generated. 

Now consider the case in which n and k are not relatively prime. We show that the assumption 
that R is deep leads to a contradiction. Thus, assuming that R is deep implies that there is a value 
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m such that R can be reaUzed as Dk,n,m = {ifn mod n : z = 0, 1, . . . , A; — 1}„. This in turn imphes 
that there exists an integer ^ such that km = in+ 1, that is, in = 1 (mod k). However, for this to 
happen, n and k must be relatively prime, a contradiction. 

Thus we have shown that R is Erdos-deep if and only if k and n are relatively prime. □ 
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